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ABSTRACT 


Calculational procedures for Monte Carlo analyses and W resonance absorption are 
presented. The widely used flat spatial neutron source assumption is shown to overesti- 
mate the effective resonance integral of the large scattering resonance in W*®® by 

9 

21 percent at a surface to mass ratio (S/M) of 0. 5 cm /g. Resonance overlap in natural 
W decreases the effective resonance integral by 18 to 7 percent over a S/M range of 

9 

0. 16 to 4. 0 cm /g relative to calculations which ignore resonance overlap. Spatial self- 
shielding in a complex geometric cell is evaluated through comparison of Monte Carlo 
calculations with exact geometric and homogenized absorber region representations. 

The code EPIGRAM, which selects energy points for resonance cross sections, is re- 
ported. 
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APPLICATIONS OF MONTE CARLO ANALYSIS TO TUNGSTEN 
RESONANCE ABSORPTION 
by Robert M. Westfall 
Lewis Research Center 

SUMMARY 

Accurate calculation of capture rates in tungsten isotopes is required for reactors 
using tungsten as a structural or shielding material. Widely used analytical procedures 
employ approximations in calculating resonance escape probabilities for arbitrary geo- 
metric configurations or for calculating effective resonance integrals. The accuracy of 
these approximations is evaluated for three cases by the Monte Carlo method. The 
first case estimates the validity of the flat spatial neutron source assumption used in 
computing collision probabilities for absorbers as applied to the large 18. 83-eV scatter- 
ing resonance in tungsten 186. The use of the flat source approximation results in an 
overestimation of resonance absorption that varies widely with absorber lump size and 
with the magnitude of scattering and absorption resonance half widths. 

The second case demonstrates the effect of ignoring resonance overlap among the 
isotopes of natural tungsten. Comparison of the calculated values with and without reso- 
nance overlap show that neglect of overlap for natural tungsten overestimates the effec- 
tive resonance integral by 18 to 7 percent over a surface to mass ratio range of 0. 16 to 
4. 0 square centimeters per gram. 

The third application is a demonstration of how the detailed geometry treatment by 
the Monte Carlo method can significantly affect the resonance absorption in a complex 
cylindrical tungsten fuel element. Monte Carlo calculations of the resonance escape 
probability for such a fuel element show a significant difference in absorption probabili- 
ties between calculations which represent the absorber materials in discrete geometry 
or as a homogenized central smeared region. 

As part of the procedures, presentation is made of the computer code EPIGRAM 
which provides zero temperature, Breit-Wigner, single level cross sections (the cross 
sections are Doppler broadened with another code) at those energy points from which the 
intermediate cross sections can be interpolated within prescribed absolute and percent- 
age deviations. The resultant minimum number of cross sections at assured interpola- 
tional accuracies improve both the speed and the reliability of subsequent Monte Carlo 
calculations using the DRAMA code. 


INTRODUCTION 


The refractory properties of tungsten (W) make it a desirable material for use in 

high temperature reactors. However, the presence of large neutron absorbing reso- 

1 82 1 88 

nances and their associated high thermal cross sections in the isotopes W , W , and 
1 88 

W 1 D precludes the use of natural tungsten in thermal reactors for nuclear propulsion. 

184 

It has been suggested (ref. 1) that tungsten enriched in the isotope W , which has a 

small resonance absorption integral and a low thermal absorption cross section, would be 

184 

applicable. Since the cross sections of a mixture highly enriched in the W isotope 
would be dependent upon the residual absorbing isotopes, an accurate method for calcu- 
lating resonance absorbtion in all isotopes is required. 

Some widely used methods of treating resonance absorption (refs. 2 to 4) have been 
applied with accuracy to materials with widely spaced, predominantly absorbing reso- 
nances such as uranium 238. Under these conditions, the analytical treatment conven- 
iently and validly assumes a flat spatial neutron source through the absorber region to 
compute collision probabilities and also complete energy flux recovery between adjacent 
resonances. These assumptions may not apply to tungsten because several of the tung- 
sten resonances are closely spaced and have large resonance scattering widths. The 
Monte Carlo method can calculate resonance absorption without requiring these assump- 
tions since this method has the advantage of following the slowing down of the neutrons 
as discrete particles. 

The procedures described in this report have required the compilation of resonance 
parameters for the tungsten isotopes, the development of computer codes for the calcu- 
lation of cross sections, and the calculation of effective resonance integrals and reso- 
nance escape probabilities with the Monte Carlo code, DRAMA (ref. 5). 

Three cases of resonance absorption are studied with the Monte Carlo method. The 
first case is an evaluation of the error associated with the use of the flat spatial neutron 

source assumption in analyzing absorption by the large 18. 83 -eV scattering resonance of 
1 88 

W . The second case is an evaluation of the effects of resonance overlap present in 
natural tungsten over a surface to mass ratio range of 0. 16 to 4. 0 square centimeters 
per gram. The third case studied is that of spatial self-shielding in a complex cylindri- 
cal tungsten fuel element. Detailed consideration of exact cell geometry relative to a 
homogized absorber region representation results in a significantly lower value of the 
cell resonance absorption probability. 

The appendix contains a description of EPIGRAM, a computer code developed to 
identify the energy points required to provide the minimum number of resonance cross 
sections at given interpolational accuracies. Additionally, the code provides the cross 
section arrays and the resonance absorption integral at infinite dilution. Operating in- 
structions, a FORTRAN listing, and a sample problem are presented. 
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TUNGSTEN CROSS SECTIONS 


w 


Resonance Parameters 

The accuracy of the Monte Carlo calculation is dependent upon the precision of the 
cross sections of the constituent materials and the completeness with which the cross 
sections are sampled over the energy spectrum. The resonance parameters for the 
tungsten isotopes have been remeasured in the last several years (refs. 6 to 10). Also, 
recent measurements on isotopically enriched tungsten samples have provided more ac- 
curate values of the isotopic -thermal absorption cross sections (ref. 11) and the infi- 
nitely dilute resonance integrals of and (ref. 12). 

The compilation of tungsten resonance parameters (listed in the natural tungsten 
sample problem in the appendix) has been drawn primarily from the Oak Ridge National 
Laboratory - Rensselaer Polytechnic Institute measurements of the last few years. The 
main source is reference 6 from which parameters for all resonances above 87. 4 eV are 
taken. Values for the low energy resonances were selected on the basis of consistency 
with reference 6. The sources of the low energy data are as follows. The first two res- 

1 ft? 1 o o 183 

onances in W and the first resonance in W are from reference 7. For W , the 
second resonance is from reference 8, the third from reference 6, and the fourth, fifth, 
and sixth are from reference 9. 

182 184 

The bound level resonances in W and W have energies which were assigned to 

make up the differences between the measured thermal cross sections of reference 11 
and the values calculated from the known parameters. The bound levels are based on the 

average neutron width and the assumed capture width of reference 6. 

186 

Parameters for the important first resonance in W are taken from reference 12 
which compared absorption integrals calculated with various capture half widths with a 

precise measurement of the dilute resonance integral. 

1 80 

W has a very small natural abundance (0. 0013); and it was not until an enriched 
sample (0. 0693) became available that experimental measurements could be made. Thus 
the five resolved resonances of reference 10 constitute all the present information on 
this isotope. 

Calculations of the thermal cross sections and dilute resonance integrals based on 
these parameters have been made for the tungsten isotopes. These are tabulated in the 
sample problem in the appendix and are compared with experimental values. 


Data Preparation 


In the Monte Carlo calculation, there is a choice of calculating Doppler broadened 
cross sections after each neutron collision or of energy interpolating the values from 
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tables of precalculated cross sections. Of the two methods, the interpolation technique 
is much faster. Since a typical calculation contains several hundred thousand neutron 
collisions, time becomes an important factor and the interpolation scheme is usually 
employed. The main area for concern with the interpolation scheme is that the tabulated 
energies must be chosen such that the cross sections can be interpolated between them 
with the desired degree of accuracy. For this purpose the computer code EPIGRAM was 
developed. 

This code, which is discussed in detail in the appendix, provides zero temperature, 
Breit-Wigner, single level cross sections. The cross sections are Doppler broadened 
through the use of the DBCS code described in reference 13. Tables of cross sections 
for the separate nuclides, along with the corresponding energy tables, are entered into 
the DRAMA code for the Monte Carlo calculation. After each neutron collision, the 
DRAMA code ascertains the interacting nuclides, interpolates the cross sections semi- 
logarithmically in energy (see the appendix), and multiplies the cross sections by the 
proper number densities. Thus, the probabilities for the various events over the neu- 
tron energy spectrum are provided by a method which has sufficient speed, desirable 
flexibility, and assured accuracy. 


CALCULATIONAL PROCEDURE OF THE DRAMA PROGRAM 

The DRAMA program uses Doppler broadened cross sections to calculate the reso- 
nance escape probabilities of configurations of interest. Since it directly affects reactor 
criticality, the resonance escape probability p is an important parameter of the fuel 
cell configuration. The resonance escape probability is a function of the effective reso- 
nance integral. For a given material, the effective resonance integral I is a meas- 
ure of the integrated absorption probability over a specified energy range. Since the ef- 
fective resonance integral varies primarily with the surface to mass ratio, it is an im- 
portant parameter with which absorbing materials can be compared. 

The procedure followed in this study has been to calculate the resonance escape 
probabilities of either complex fuel cell configurations or simple absorber geometries, 
obtaining the effective resonance integral through the following expression: 

( ~ N A V A l ett\ 

p = exp I — A _ eH (1) 

\ V M^ S s / 

where N A is the atom density of the absorber, £2 g is the average slowing down power 
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of the moderator, and and are the volumes of the absorber and moderator, 

respectively. 

Calculational Geometry 

The DRAMA program considers a geometry which is based upon an absorber- 
moderator cell. Upon intersecting the outer cell boundary, the neutrons T coordinates 
are adjusted so that they re-enter the cell in the same relative position that they would 
enter the adjacent cell of an infinite repetitive array. The outer cell boundaries may be 
cylindrical, rectangular, or hexagonal. The cells may be bounded by planes in the ver- 
tical direction. The internal cell boundaries, separating regions of constant material 
composition, can take any shape described by a first or second degree equation. The 
present version of the code can consider up to eight materials in any of the 30 internal 
regions. The program is limited to 15 materials overall, six of which can be resonance 
absorbers. A further limitation is the total of 50 internal and external boundaries used 
to describe the cell. 

Resonance escape probabilities have been calculated for complex two-dimensional 
hexagonal cells containing concentric cylindrical regions, an example of which is shown 

Uranium 238 

Enriched uranium- 
aluminum alloy 


Aluminum support 
mandrels surrounded 
by natural tungsten foil 

tungsten 

Figure 1. - Cross-sectional view of complex cylindrical tungsten fuel element. 
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in figure 1. Radial and axial leakage can be calculated by placing thick regions of 
"black” absorbing material adjacent to the outer cell boundaries. 

When the calculational objective was obtaining effective resonance integrals, the 
slab lattice model shown in figure 2 was used. This simple model is easy to set up and 



Half cell for 
calculation 


Figure 2. - Slab cell model in 
semiinfinite lattice. 


does not require long running times for results of acceptable statistical quality. Sym- 
metry conditions allow calculation of the half cell. Reflective cell boundary conditions 
restrict the surface of the absorber to a single plane. For a given absorber density p a 
( g/cm\ the absorber thickness r (cm) is related to the desired surface to mass ratio 

O 

S/M (cm /g) through the following expression: 


T = 


2 

P A ( S/M) 


(2) 


A procedure used to reduce running time was to minimize the moderator thickness. 
Minimum moderator thicknesses for maintaining an unperturbed neutron flux distribu- 
tion, that is, a 1/E neutron flux, were established for several absorber material con- 
figurations. 


Neutron Slowing Down Treatment 

The DRAMA program employs a rejection technique to select the initial neutrons 
from an assumed l/E neutron flux distribution over the energy interval between the 
maximum energy for the calculation and the minimum energy resulting from a single 
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elastic collision with the moderator atom. 

Since inelastic scattering and p-wave and higher l - wave resonances are ignored in 
the calculation of the cross sections, the maximum energy should not exceed several 
thousand electron volts. An initial statistical weight of one is assigned to each neutron. 
This is reduced by the absorption probability over the neutron path length as it is scat- 
tered down in energy. The scattering is isotropic in the center of mass coordinate sys- 
tem. Each history is followed until its statistical weight falls below 0. 00001 or the neu- 
tron’s energy falls below the cutoff for the calculation. Since upscattering is ignored, 
the cutoff energy should be in excess of 5 kT where k is Boltzmann's constant and T 
is the temperature of the medium in °K. 


Statistical Accuracy 

Calculation of the probable error has been used to determine the statistical accuracy 
as a function of the number of neutron histories. The resonance escape for the Monte 
Carlo calculation is defined as the sum of the terminal neutron weights clk divided by 
the number of neutron histories H. The variance of the resonance escape V(p) is the 
second moment of p minus the square of the first moment; that is, 



i 


The probable error PE(p) is calculated as, in reference 14, 


PE(p) = 0. 6745 \l— M. 

jr H - 1 


(3) 


(4) 


There is a 50 percent expectation that the correct resonance escape value lies within 
±PE(p). To increase the expectation, the appropriate coefficient taken from a table of 
the normal distribution function must be substituted. For ten thousand neutron histories, 
the probable errors of the resonance escape values run from 0. 5 percent for high p 
values (p > 0. 95) to 1 percent for low p values (p < 0. 85). Since the effective reso- 
nance integral has a logarithmic dependence upon the resonance escape, the associated 
probable errors vary in the opposite direction. Typically, they vary from 2 percent for 
thin, moderately absorbing samples to 1 percent for thick, highly absorbing samples. It 
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should be noted that the thick samples have lower effective resonance integrals so that in 
absolute terms their associated probable errors are usually much less than those of thin 
samples. This indicates the increased difficulty of adequately sampling thin absorbers. 


THE FLAT SPATIAL NEUTRON SOURCE APPROXIMATION 

The flat spatial neutron source approximation, as applied in the Nordheim integral 
method of reference 2 to calculate collision probabilities, was early recognized to cause 
an overestimation of the resonance absorption for thick lumps of material with high ra- 
tios of resonance scattering to absorption (ref. 15). The validity of the flat source ap- 

IOC 

proximation for the large 18. 83-eV scattering resonance in W has been studied by 
Cohen (ref. 16) using transport calculations. 

The present Monte Carlo calculations were applied to the same problem to deter- 
mine how the effective resonance integral varies with both absorber lump size and the 
ratio of scattering to absorption. The procedure involved a comparison of ZUT (ref. 4), 
which uses the flat source approximation, and DRAMA calculations of the effective reso- 
nance integral arising from the single 18. 83 eV resonance. Consideration of the single 
resonance precludes the effect of resonance overlap as a complicating factor. A range 



Figure 3. - Scattering cross section against energy for 
18. 83 eV resonance. Capture reaction half width 
r = 0.052 eV ; temperature T = 300° K. 
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Figure 4. - Absorption cross section against energy for 
18. 83 eV resonance. Capture reaction half width 
T y = 0. 52 eV ; temperature T * 300° K. 


of S/M values from 0. 12 5 to 1. 12 5 square centimeters per gram was considered. The 
ratio of scattering to absorption was varied by using scattering half widths of 1. 0, 0. 315, 
and 0. 1 eV with a constant capture half width of 0. 052 eV. 

Both the Monte Carlo and ZUT calculations were based on the slab cell geometry. 

In the Monte Carlo calculation, 4 inches (10. 16 cm) of water separating the absorber re- 
gions maintained a 1/E flux incident on the slabs. The atom density of W 1 D in the 
slabs was based upon a physical density of W ° of 19. 3 grams per cubic centimeter. 

The energy range considered was between 100 and 8. 07 eV (the 3/7 E Q cutoff of the ZUT 
calculation). 

The 600 energy point, Breit-Wigner, single level cross sections of the Monte Carlo 
calculation were Doppler broadened to 300° K with the same shape functions used in the 
ZUT calculation. The scattering and capture cross sections for the actual and for ex- 
treme neutron widths are plotted in figures 3 and 4. It is seen that for the scattering 
cross section, both the peak value and the resonance width vary directly as the neutron 
width. Although the peak value of the capture cross section varies inversely as the neu- 
tron width, the integrated absorption probability over the resonance varies directly. 
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TABLE I. - COMPUTED RESONANCE INTEGRALS OF W 186 
BETWEEN 100 AND 8. 07 eV 


[Resonance energy E q = 18. 83 eV; capture reaction half width 
T = 0. 052 eV; temperature T = 300° K; density 
p = 19.3 g/cm 3 .] 


Scattering 

Type of 

Surface to mass ratio, S/M, 

reaction 

calculation 


crcr/g 


half width, 

r n> 

eV 


0. 125 

0. 5 

1. 125 

0. 1 

I (ZUT) 1 ” 

9.83 

20. 63 

30. 52 


I (Monte Carlo) ^ 

9.05±0.35 

19. 68±0. 98 

30. 02±1. 8 


i /i * 
z' me 

1.09 

1. 04 

1.02 

0.315 

I (ZUt/ 

15.58 

33. 96 

51. 08 


I (Monte Carlo/ 

12. 56±0. 43 

28. 02±1. 14 

45. 8±2. 2 


I /I * 
z' me 

1.24 

1.21 

1. 11 

1.0 

i (zut/ 

21.6 

52. 98 

81. 71 


I (Monte Carlo/ 

12. 97±0. 45 

42. 45±1. 55 

64. 34±3. 1 


1 /I * 
z' me 

1. 67 

1.24 

1.27 


^Resonance integral in barns. 
\A mc = I (ZUT)/I (Monte Carlo). 


The results of the study of the flat neutron source approximation are given in table I. 
Each of the Monte Carlo values represents 10 000 neutron histories, and the probable 
errors shown reflect their statistical accuracy. For a neutron width of 0. 1 eV (on the 
same order as the capture width), the ZUT values do not seriously overestimate the ef- 
fective resonance absorption integral. However, the dependence of the overestimation 
on absorber lump size is clearly seen. At the neutron width of 0. 315 eV which corre- 
sponds to measured values of the actual resonance, the overestimation is seen to be in 
excess of 20 percent in the range of S/M values which are found in many reactor appli- 
cations. The extreme case of the neutron width of 1. 0 eV illustrates very large over- 
estimations of I e ££ by the ZUT calculation. 

Cohen in reference 16 compared transport calculations employing a fine spatial 

mesh and 48 energy groups with calculations done with the ZUT code. Analysis was 

1 

done on the effective resonance absorption integral of the same 18. 83-eV level in W 
using a scattering half width of 0. 317 eV and an energy range of 10. 38 to 27. 85 eV. An 
effective resonance absorption integral overestimation of 13 percent at S/M = 1. 0 
square centimeter per gram agrees closely with the 11 percent value at a S/M value 
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of 1. 12 5 square centimeter per gram in table I. However, with the substitution of a hy- 
perbolic cosine neutron distribution in the absorber region to compute the escape proba- 
bility values as used by the ZUT code, Cohen has shown that the Nordheim integral 
method can be used to accurately calculate resonance absorption in thick slab geometries 
of highly scattering materials. 


RESONANCE OVERLAP IN NATURAL TUNGSTEN 

A commonly used assumption is that resonances are sufficiently spaced in energy to 
allow the complete recovery of the incident neutron flux to a 1/E distribution and thus 
do not interfere with each other. Natural tungsten has several closely spaced, large 
resonances in which overlap does occur. The present Monte Carlo calculations were 
made to determine the magnitude of the overlap effect in natural tungsten for a series of 
Vs/M values from 0. 4 to 2 centimeters per gram*^ (S/M varies from 0. 16 to 
4. 0 cmVg). 

The effective resonance integrals were computed for samples of natural tungsten and 
for the four principal tungsten isotopes, W^, and at natural abun- 

dances. The overlap effect was estimated from the comparison of the value calculated 
for the natural tungsten with the sum of the values calculated for the constituent isotopes 
separately. 

A repetitive slab lattice of alternate tungsten- water regions was used in the calcula- 
tions. Two inches (5. 08 cm) of water was sufficient to eliminate cell interaction effects. 
Six hundred energy point sets of Breit-Wigner, Doppler broadened (300° K) cross sec- 
tions described each isotope. The energy range considered was from 2200 to 0. 5 eV. 

The cross sections were calculated using resonance parameters available prior to the 
publication of reference 6; however, these are essentially the same as those found in the 
appendix. Atom densities for the isotopes were based on a tungsten density of 
19.3 grams per cubic centimeter. Each calculation included 10 000 neutron histories 
except for the thin sample case, Vs/M = 2. 0 centimeters per gram^, which required 
20 000 histories for acceptable statistical accu racy. The probable errors of the Monte 
Carlo met hod v ary from 1. 3 percent for Vs/M = 0. 4 centimeter per gram*^ to 2. 5 per- 
cent for Vs/M = 2.0 centimeters per gram^. 

Figure 5 shows the Monte Carlo effective resonance integrals for natural tungsten 
as calculated for the natural isotopic mixture and as the sum of contributions of sepa- 
rated, noninteracting isotopes. Also shown in figure 5 are experimental epicadmium ef- 
fective resonance integrals (ref. 17). The Monte Carlo natural mixture results and the 
experimental results are seen to be in reasonable agreement. Inclusion of self-shielded 
unresolved resonance contributions of 0. 5 to 1. 5 barns as calculated by the TUZ code 


11 



Absorber lump dimension, Vs/M, Vcm^/g 

Figure 5. - Natural tungsten effective resonance 
integrals. 

(ref. 4) in the Monte Carlo values for the range of VS/M improves the agreement 
slightly. 

Comparison of t he Mo nte Carlo values in figure 5 shows the overlap e ffect to vary 
from 18 percent at VS/M = 0. 4 centimeter per gram^ to 7 percent at Vs/M = 2. 0 
centimeters per gram^. By tabulating the absorptions between prescribed energy 
points in the Monte Carlo ca lculat ions it was possible to determine where the overlap 
occurs. For all values of Vs/M, 80 percent or more of the overlap effect is attribut- 

IOC 1 op 

able to the interaction of W and W which have large resonances at 18. 83 and 
21. 09 eV, respectively. For the thicker samples, less than 10 percent of the overlap 
effect occurs in the 1/v region above 0. 5 eV. 

The Nordheim integral method has been extended in the GAROL code (ref. 18) to ac- 
count for resonance overlap in calculating effective absorption cross sections. Calcula- 
tions of W^, and between 2. 38 and 61. 4 eV show overlap effects of 6. 5, 

1. 2, and 9. 3 percent, respectively (ref. 19). 

SPATIAL SELF-SHIELDING OF RESONANCE ABSORPTION IN A TUNGSTEN 
FUEL ELEMENT OF COMPLEX CYLINDRICAL GEOMETRY 

This section illustrates the calculation of the resonance escape probability for a 
complex cylindrical fuel element and demonstrates that the spatial self-shielding can be 
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precisely accounted for by the Monte Carlo method. The capability of the DRAMA code 
to represent the internal and external boundaries of the lattice cell has already been dis- 
cussed. This capability, coupled with the Monte Carlo method’s continuous treatment of 
the phase space, permits an exact geometrical representation in the neutron slowing 
down calculation. Thus the Monte Carlo method has an inherent advantage over more 
approximate methods in calculating the spatial self- shielding of resonance absorption in 
a multiple region cell. The present study is of an experimental fuel element (ref. 19). 
The effects of the spatial self-shielding in shown through comparison of Monte Carlo 
calculations containing discrete and smeared representations of the annular absorber re- 
gions. The resonance escape value derived from a multigroup neutron transport calcula- 
tion of the same fuel element is discussed qualitatively in the light of the Monte Carlo 
findings. 

The fuel element geometry is shown in figure 1. Detailed compositions and dimen- 
sions of the material regions can be found in reference 19. Briefly, the overall pitch of 
the lattice cell is 7. 62 centimeters. The 0. 16 5- centimeter outer aluminum pressure 
tube has an outer diameter of 6. 5 centimeters. Resting on five aluminum support man- 
drels and an inner support tube are six subassemblies. The inner five subassemblies 

have 0. 0127 centimeter of natural tungsten surrounded by 0. 1067 centimeter of 35 weight 

235 

percent uranium (93. 5 percent U ) - aluminum alloy. The outermost subassembly has 
0. 762 centimeter of natural tungsten surrounded by 0. 1016 centimeter of depleted ura- 
nium (0.22 percent U^^). 

The geometrical model used in the Monte Carlo calculations was the same as that 
shown in figure 1 with the exception that the inner void was smeared into the thick alu- 
minum support tube. This results in a 30 region geometry, the present capacity of 
DRAMA. Calculations were of the quarter cell in the x-y plane. Each of the six reso- 
nance absorbers (four isotopes of tungsten, two of uranium) were entered discretely. 

They were represented by 600 energy point tables of Breit-Wigner, single level cross 
sections Doppler broadened to 300° K. The calculations included a slowing-down energy 
range from 2034. 68 to 0. 532 eV. 

The resonance escape probability for the discrete geometry model was calculated to 
be 0. 708. This value represents 6000 neutron histories and it has an associated prob- 
able error of 0. 003. The smeared geometry contained two regions, the water moderator 
and the remainder of the cell. The resonance escape probability after 6000 histories 
was 0. 691 -with a probable error of 0. 002. The absorption in the cell during slowing 
down is therefore 5. 5 percent greater for the smeared case than for the case with the 
discretely represented geometry. 

Current methods of reactor analysis involve the use of the multigroup S n transport 
calculation to determine the reactor multiplication factor. The general procedure is to 
consider the fuel element and its surrounding moderator as a cell in an infinite lattice. 
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Material cross sections are prepared with cross section averaging programs such as 
GAM II (ref. 20). Using an S n transport program such as TDSN (ref. 21) and the mate- 
rial cross sections, a cell calculation is performed to obtain flux weighted cell cross 
sections. The flux weighted cell cross sections are then used in a transport calculation 
of the full reactor core to determine the multiplication factor. 

The cell transport calculation corresponds to the Monte Carlo calculation described 
previously. It is of interest to derive the resonance escape probability from the cell 
calculation and compare it to the Monte Carlo value. However, only a qualitative com- 
parison can be made because of several approximations inherent in both the cross sec- 
tion averaging and the cell calculation. 

A 16 energy group S^, Pj, one-dimensional transport calculation was done for the 
same cell geometry as that used in the discrete region Monte Carlo calculation. Ninety- 
nine mesh points were distributed over the 30 discrete material regions. The hexagonal 
outer cell boundary was approximated by the circle which conserves the area of the mod- 
erator. 

The 16 energy groups of the material cross sections included four groups in the res- 
onance energy region of the Monte Carlo calculations. The GAM II program was used to 
calculate 34 fine group resonance cross sections and flux weight them into the four broad 
■groups. The calculation of the fine group cross sections is done with the ZUT-TUZ pro- 
cedure. Three approximations are involved in this initial stage of the cross section 
preparation. The absorber materials are treated as a lumped region and thus spatial 
self-shielding is ignored. The escape probability table for the lumped absorber is based 
on the flat spatial neutron distribution approximation for a solid cylindrical rod. And 
each resonance nuclide is treated separately with resonance overlap not taken into ac- 
count. The 34 fine group cross sections were flux weighted over a spectrum resulting 
from the interaction of fission spectrum neutrons with the cell materials to obtain the 
four broad group cross sections. Spatial self-shielding remains to be treated in the cell 
calculation. 

The S n cell calculation, in principle, could exactly account for spatial self- 
shielding if it contained an adequate number of mesh points. How well the 99 mesh 
points account for the spatial self-shielding cannot be determined quantitatively because 
of the approximations in the cross section preparation discussed previously and also be- 
cause of several differences between the transport and Monte Carlo calculations. 

The significance of the resonance overlap effect was shown by making Monte Carlo 
calculations of each resonance nuclide separately in the discrete region geometry. The 
resulting value of the resonance escape probability 0. 608, reflects an absorption prob- 
ability which is 39. 2 percent greater than the value with resonance overlap taken into 
account. 

The value of the resonance escape probability derived from the transport calculation 
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is 0. 669. The difference between this value and the discrete region Monte Carlo value of 
0. 708 cannot be attributed entirely to resonance overlap and spatial self-shielding. In 
addition to the several effects discussed previously, minor effects could arise from the 
circular outer cell boundary approximation or the difference in the source neutron flux 
slowing down into the resonance region. 

The Monte Carlo calculations have shown that the detailed representation of geom- 
etry results in an absorption probability which is 5. 5 percent less than that of the 
smeared cell. It is probable that the transport cell calculation does not account for all 
of this effect. Significant differences have been shown between the S n and Monte Carlo 
resonance absorption probabilities. The number of approximations in the S n transport 
calculation demonstrates the desirability of the exact Monte Carlo method for cells of 
complex geometry. 


CONCLUSIONS 

The Monte Carlo method as applied to resonance analysis can account directly for 
effects which other methods treat approximately. The results of the Monte Carlo studies 
of three of these effects upon tungsten resonance absorption are summarized here. 

The analytic approximation based on the assumption of a flat spatial neutron source 
distribution through the absorber leads to overestimations of resonance absorption in 
thick lumps of high resonance scattering material. The effective resonance absorption 

1 ftfi 

integral contribution of the large scattering resonance at 18. 83 eV in W is overesti- 
mated by 21 percent at a surface to mass ratio of 0. 5 square centimeter per gram when 
the flat spatial neutron distribution is employed. 

Resonance overlap in natural tungsten reduces the effective resonance integral at 
low surface to mass ratios relative to calculations which superpose the isotopic capture 
rates. The overlap effect varies from 18 to 7 percent over a surface to mass ratio 
range of 0. 16 to 4. 0 square centimeters per gram. 

Spatial self -shielding decreases the resonance absorption probability of a complex 
cylindrical tungsten fuel element by 5. 5 percent. The resonance escape probability de- 
rived from a neutron transport calculation of the same fuel element is significantly less 
than the Monte Carlo value. The difference is attributed to a number of effects of which 
resonance overlap and spatial self -shielding are probably the most important, but the ef- 
fects cannot be separated in the transport calculation. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, February 1, 1968, 

129-02-04-03-22. 
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APPENDIX - COMPUTER CODE EPIGRAM 


General Description 

The purpose of the code is to generate zero temperature, Breit-Wigner, single level 
cross sections at those energy points which allow interpolation of intermediate cross 
sections at a desired degree of accuracy. The minimum number of energy points thus 
produced reduces computation time. 

The procedure followed by the code is to cover the resonances with a large number 
of closely spaced energies. The corresponding cross sections are calculated. Then, by 
using the same interpolation technique applied in the Monte Carlo calculation, those en- 
ergies at which the cross sections can be interpolated within the desired margins of ab- 
solute and percentage accuracies are eliminated. 

The initial energy grid is generated upon the resonance base energies, the energies 
at one-half and two total widths off the base, and the wing points at which the total cross 
section is less than 0. 0005 of the maximum. If the wing points of adjacent resonances 
overlap, they are reset at equal intervals between the corresponding points at two total 
widths off the base energies. The initial grid is expanded by equally spacing 19 addi- 
tional points between each of the original points. The process results in a grid of 
140 N + 21 energy points where N is the number of resonances. 

Through iterative procedures, the code can consider a material containing as many 
as 1200 resonances. Ten is the maximum number of isotopes. Each of the two statisti- 
cal weight factors g associated with each isotope can have as many as 60 resonances. 

If it is desired to have more than 60 resonances associated with one g factor, the iso- 
tope can be entered twice with the second entry containing the additional resonances in 
one or both of its g factors. In this instance, care must be taken to ensure the correct 
potential scattering contribution for the isotope by using a zero value for the potential 
scattering multiplying term PSTM of the second entry. 

The code executes the interpolation and elimination-through- comparison scheme on 
the table of total cross sections. The zero temperature, Breit-Wigner, single level 
cross sections with interference between resonance and potential scattering taken into 
account are calculated according to the following formulas: 


16 



NI 

^ 

in=l 


4ttB. 


m 


A in + 1. 008986 J 



NRk 

EE 

L k=1 

i=i 

! 




J °, 'J 




^ E e _i5 sin 6 


> (Al) 


NI 2 

E /a. + 1.008986\ 2 \ N \ ' r ,r A_ , 

" B in H 5 — r ) %*e ) ^ ’’ , 

V in / < E - E o j> - (r / 2) ‘ 

‘ ' k=l j=l 


in=l 


NI 


(A2) 


NR 


k 


a f (E) = 


E /A. + 1. 008986\ 2 \ N \ S 

L 6 " E L 

in=l k=l i=l 


F n, j r f, j^r, j 


L_j (e - e 0> / + (r/2 ) 2 

i=i 


(A3) 


where 

A atomic mass 

B fractional atomic abundance 

E neutron energy, eV 

E q resonance energy, eV 

g statistical weight factor 

i yTi 

NI number of isotopes 

NR number of resonances 
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i 



cm 


R 

r 

6 

a. 


nuclear radius, 1. 5x10" 13 A 1 / 3 , 
reaction half width, eV 
-2. 196979X10 9 R Vi 

2 - 86x10 ' 9 , cm 

2? r Vi 


a cross section, b 

and the subscripts are 
c, y capture 

E neutron energy applies 

f fission 


in isotope number 

j resonance number 

k g factor number 

n, s scattering 

r resonance energy applies 

The interpolation may be done either semilogarithmically or linearly in energy. The 
semilogarithmic scheme requires more time but results in about 25 percent fewer points 
for the same degree of accuracy. The shorter table lookup time in the Monte Carlo cal- 
culation more than makes up for the longer interpolation time. Whichever scheme is 
selected, assured accuracy of the cross sections depends upon the consistency of their 
generation and application. 

For energy points i, j, and k, the cross section at point j is interpolated as 
follows: 

S emilogarithmi c: 


Linear: 



/a.\( E k- E 3)/( E k" E i) 

V a J 




(Ej - Ej) + 


CT i 


(A4) 
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Provision is made to eliminate progressive series of points j by advancing the index k 
and reinterpolating and comparing all points j. When a point j fails the tests, the k 
point corresponding to the preceding successful series of interpolations is retained to be 
printed out in the final array. Then the i index is advanced to this k value and the 
process continues. 

When a point fails on its initial interpolation, that is, it must be retained, additional 
tests are made on the percent deviation to determine if the initial grid is too coarse be- 
tween points i and k. The upper limit or the percent deviation is set at 10 percent ex- 
cept at those points where Uj + < 2 barns. At these points deviations up to 

20 percent are allowed. If these upper limits are exceeded, intermediate points at suc- 
cessively finer intervals are inserted between points i and k. Cross sections are cal- 
culated and the interpolation tests are run until all of the intermediate points fall within 
the upper limit criteria. These points are retained, the i index is advanced to k, and 
the elimination of points from the original lattice is resumed. 

Additionally, certain points of particular interest are retained. These are the max- 
imum, minimum, and inflection points associated with the resonances. Inflection points 
are detected by monitoring the sign of the interpolated deviation. Maxima and minima 
are detected by comparing successive cross section values. Neither criteria is applied 
if a point is retained for exceeding the interpolated deviation. 

Flexibility has been provided in the code to permit three irregular situations. The 
first is the inclusion of the contributions from bound level resonances at negative reso- 
nance energies. These are included in the input in the normal manner with the number 
of negative resonances NNR specified on the first control card. This results in their ex- 
clusion from the initial energy lattice, however, the subsequent cross section calcula- 
tions do account for their contributions (by taking the absolute values of the resonance 
energies). 

The second situation is one in which resonances of interest lie above or below the 
energy range of the problem. Resonances below the low energy cutoff EMN can be ac- 
counted for by using the NNR variable in the same manner described previously. Reso- 
nances above the high energy cutoff can be treated in the same way by specifying their 
total number NRSAEX on the first control card. 

The third situation arises from resonances for which the spin state has not been de- 
termined. An average g factor value of one-half is used in resolving their parameters. 
This procedure has led to three effective g factors for these isotopes. Using the same 
technique as that employed in including more than 60 resonances for a single g factor 
value will allow the inclusion of these resonances; that is, by entering the isotope twice, 
using the g = 1/2 resonances in the second entry, and setting its potential scattering 
multiplying term equal to zero. An application of this procedure appears in the sample 
case. Use of the g = 1/2 resonances represents an approximation, and it can result in 



excessive interference and small negative scattering cross sections in the s-wave dip of 
small resonances. Provision has been made to set these cross sections equal to 
0. 01 barn. The procedure is more accurate than ignoring the resonances or arbitrarily 
assigning spin states. 

Finally, the code calculates the incremental contribution to the resonance absorption 
integral at infinite dilution. The contribution between energies Ej and E 2 is calcu- 
lated by applying the linear interpolation expression for the absorption cross section to 
the following formula: 



The printed output includes the retained energies with their cross sections and the 
corresponding dilute resonance integral. These are followed by an edit which contains 
information on how many of the points are retained under the various criteria and the 
average percent deviation of the discarded points. The energies and their cross sections 
are punched out on cards. 

The program is written in the FORTRAN IV language. Running times on an IBM 
7094 n vary from 0. 1 minute for a five resonance problem to 15 minutes for a full en- 
ergy grid containing 115 resonances. If a number of the resonances are above or below 
the energy grid, such as in the sample case, the running time is considerably reduced. 

A 27 resonance problem takes less than 2 minutes. 


Input Instructions 


I Problem Identification Cards 

Alphanumeric information in card columns 3 to 72 will be printed previous to prob- 
lem output. The last card must have a number in card column 1. 
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II Control Card 


Quantity 

Format 

Card columns 

Remarks 

N 

no 

1-10 

Total number of resonances 

NI 

no 

11-20 

Number of nuclides 

FD 

F 10. 5 

21-30 

Fixed difference between calculated and interpo- 
lated values allowed 

PD 

F10. 5 

31-40 

Percentage difference between calculated and in- 
terpolated values allowed 

EMN 

F10. 5 

41-50 

Minimum energy for output grid 

EMX 

F 10. 5 

51-60 

Maximum energy for output grid 

INTYP 

15 

61-65 

= 0 Semilogarithmic interpolation will be used 
= 1 Linear interpolation will be used 

NNR 

15 

66-70 

Number of resonances below EMN 

NRSAEX 

15 

71-75 

Number of resonances above EMX 

PCH 

F5. 0 

76-80 

= 0. 0 No punched cards 


= 1. 0 One set per card of energy, total, scatter- 
ing, capture, and fission cross sections in 
Format 1P5E 12. 5 


III Data Card For Each Isotope (NI cards) 


Quantity 

Format 

Card columns 

Remarks 

A 

F 10. 5 

1-10 

Atomic mass number 

AB 

F 10. 5 

11-20 

Fractional abundance 

R 

Ell. 5 

21-31 

Nuclear radius (cm) 

NGI 

19 

32-40 

Number of g factors associated with isotope 
(1 or 2) 

PSTM 

F10. 5 

41-50 

Potential scattering term multiplier (1. 0 unless 
isotope appears more than once, 0. 0 for sec- 
ond entry) 
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ii in mi m«i 


IV Data Card For Each g Factor (two NI cards) 


An extra blank card must be inserted for each isotope having only one g factor (see 
sample problem). 


Quantity 

Format 

Card columns 

Remarks 

GL 

F10. 5 

1-10 

Value of g factor 

NR 

110 

11-20 

Number of resonances associated with each g 
factor 


V Data Card For Each Resonance (N cards) 


Quantity 

Format 

Card columns 

Remarks 

EB 

F10. 5 

1-10 

Energy of resonance (eV) 

NGTG 

110 

11-20 

Identification of g factor of resonance, g fac- 
tors numbered in order of entry 

GMNP 

F10. 5 

21-30 

Neutron scattering half width (eV) 

GMGP 

F10. 5 

31-40 

Radiative capture half width (eV) 

GMFP 

F10. 5 

41-50 

Fission half width (eV) 


FORTRAN Listing 


c 

c **** Ep I GRAM **** 

C 

C ENERGY POINT IDENTIFICATION AND CRCSS SECTION GENERATION FOR 

C RESONANCE AdSURBING MATERIALS 

C 

C COMPLEX ALGEBRA PROGRAMED BY MARY ANN ARNOLD, LEtolS RESEARCH CENTER 

C 

DIMENSION ER(SCE), $GSR(5C5), SGAR(505), SGFR(505), SGTR(505), 

1ES( 275C), SG S ( 27 50) r SGASI2750), SGFS(2750>, A ( 10 ) , AB(IO), R(I0), 
2G(IC,2J, ECH2C, EC), GMN<20, 60 , GMG (2C,6C) , GMF<20 , 60 ), NR(20), 

2GL ( 20 ) f NR G ( 20 ) , GI110), COEF (10) , F(501), ET (260 ) , ST(260), 

4SAT(2£C), SF T ( 26C ) , EB(12CCi, GNNPU200), GMGP(1200), GMFP(1200), 
5N GTG( 1 2CC ) , NGK1C), ADRA I ( 505 ) , PST Ml 0 ) , TITLEQ4) 

EQUIVALENCE ( SGA S ( 1 ) , NG TG ( 1 ) ) » ( SG A S (1201 ) , G N F P ( 1 ) ) , (SGA$(2401), 

1GMGPC 1 ) ), ( SAT( 1 ) , GMGP (351) ),'(SFTH) , GMGP(611> ) 

COMPLEX DEN, EX1, Tl, T2 , CEXP1, V A 13 
WR I TE ( E, ICC ) 

C 

C REA C AND V»R I TE TITLE CARDS 

C 

74 REA C( 5,112) NTITLE, (TIT LE Cl) ,1=1,14) 

WRITE( 6,114 ) (TITLE(I), 1=1,14) 

IF(NTITLE) 74,74,72 
C 

C REAL GENERAL INPLT 
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m 


C 

7 5 REAC (MCI) M ,M ,FD ,PD ,EMN ,EMX ,1 MVP , NNR,NRSAEX ,PCH 

IF( IN T VP >76, 79,78 
76 WR I T E ( 6, 1 J5 ) 

GO TO 6C 
79 WRITEt {,116} 

C 

C REAC NUCL IDE DATA 

C 

8 C REAC ( f > 1 C 2 ) (A(I), ABU), R(I), NGI(I), PST M { I ) , 1=1, Ml 

I GN = N 1*2 
C 

C REAC G FACTOR DATA 

C 

REAC IS1C2) (GL( I ) ,NR(I ) ,1 =1 ,IGM 
C 

C INITIALIZE PARAMETERS 

C 

KOI = C 
CKO 2 = C.C 
NK ST = C 
KC = C 
NK = 1 
KZ = C 

ICU = C 
NJ = C 
C 

c 

C ASSOCIATE G FACTORS WITH NUCLIDES 

C 

C 

DO 1 I = 1 , N I 
M = 2* I 
MM =M— 1 

G ( I, 1)= GL(MM) 

1 G( I, 2 ) = GUM ) 

C 

C REAC RESONANCE DATA 

C 

REAC d,lC4) (EB(I).NGTGU ) , G MNP (I ) , G MG F ( I ) , GMF P (I ) , 1=1, N) 
C 

C ZERC COUNTER FIELD 

C 

DO 2 I = ), IGN 

2 NRG( I 1 = C 
C 

C 

C MATCH RESONANCE DATA WITH G FACTORS 

C 

C 

DO 2 I = 1 , N 
N T 1 = N GTG ( I ) 

NRC-(NTJ) =NRG(NT1) + 1 
N T 2= NRGINTl ) 

EO ( NT 1,NT 2 ) = EB( I) 

GMN ( N T ] ,N T 2 ) = GMNP ( I ) 

GMG(NT1,NT2) = GMGP ( I ) 

GM F ( N T 1,N T 2 ) = GMF P ( I ) 

3 GMNP ( I)= GMNP ( I ) + GMGP III + GMFP(I) 

C 

C ZERO INITIAL TOTAL CROSS SECTION FIELD 

C 

DO 21 1*1 , 5 7 SC 
2 1 SGASi I )=0 .C 
C 

c 

C CALCULATE NlCLIDE CROSS SECTION COEFFICIENTS 

C 

C 

DO A J =1,NI 

AF = ( ( A ( J ) + l.C08986)/A( J))**2 
A COEF(J) = AB ( J )* AF* 1 • 25 663 704E2 5 


27 

26 


39 


65 


23 


r> n o o n nor^or>or>o r> o n onn 


WRITE INPUT DATA BY NUCLIDE ,G F AC T C R f RE S CN A NC E 
DO 55 1= 1 » N I 

WRITER, 11C) 126 

AP = A( I ) 

A BP =AB( I) 

RP =R{ I ) 

ING=NG H I ) 

CO 55 JG= 1, ING 
GP = G( I, IG ) 

NF = 1*2— niG 
NGR = NR!NF) 

CO 55 IN = 1 t NGR 
EP = EO(NF, IN ) 

GMN W=GMN ( NF, IN ) 

GMGW = GMG( NF, IN ) 

GMFW = C-MFtNF, IN ) 

55 WRITE! £,111) AP, ABP, RP, GP, EF, GRNW, GMGW , GMFW 
ZERC ITERMED IATE CAPTURE AND FISSICN FIELDS 

144 

CO 25 1=1, 26C 

SAT! I )=C.C 
25 S FT ( I ) =C • C 


*** THROUGH 5, SET UP P ARRAY ON EO, i/2 + 2 GAMMA TCTAL CFF EO, 
ANC WING POINTS *** 


P! 1 ) = EMN 
ND = HNNR 
IC= 1 
J = 1 


LOOP THROUGH STATEMENTS 14,18 ON 83 RESCNANCES AT A PASS 


NOES = N - NRSAEX 
IF! ND E S «L T *83) GO TO 15 
MD = 8 2 
C-0 TO 1 6 

15 MC = NCES 

16 IC = IC + 1 

I F! ND • EQ • 1+NNR ) GO TO 17 
PI 1) = P( J + l) 

J = 1 
NSK = C 

17 DO 5 I =ND , MD 

IF! NSK .EQ . II GO TO 82 
IF! EB( D.NE.EB! I + D) GO TO 81 
IF! GMNP! I l-GMNP! 1 + 1) ) 82,83,83 

82 NSK - C 
J = J- 7 
GO TO 5 
82 NSK = ] 

8 1 GT = GMNP! I ) 

GT 2 = G T*G T 
GM = 5 • 5 
6 WS = GM*C-T 
tw = EB ( IJ-bS 
IF! EW-EMN ) 5 7, 5 7,566 
5 7 P( J + l) = EMN + C.C5 
GO TO 58 

566 CR = ( 5QR T ( E B ( I ) ) * G T 2 ) /( SQRT ( E W) * (4.0*WS**2 + GT2 ) ) 2C2 

I F! CR .LT.C.0CC5) GO TO 7 
GM = GM*5.0 
GO TO 6 


2C3 
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on nonnnnn on non o oonoooo nnooo 


7 

56 


PI J 4l ) = EB( I) - US 
P I J 42 ) = EB ( I )— 2 • C*G T 


P( J + 3 ) = EB( I )-C. 5*GT 
P(J44) = EB I I ) 

P( J+5 ) = EB C I ) + C.S*GT 
PI J46 ) = EB ( I ) + 2.C*GT 
PI J 47 ) = EB C I ) + V»S 
J = J *1 
J = J — 1 

IFIMD.NE.NDES) GU TO 6 


I FI PI J 4 1 ) .LE *E M X ) GO TO 73 
P( J + ll = EMX - l.C 
PI J +2 ) a EMX 
J = J 4 1 


ChECK AND RESET IF UI NG POINTS OVERLAP 


6 DO 9 I L = 2 » J > 

L a IL - 1 

I F I P I IL).GT.P(L)) GO TO S 
K = IL - 2 
K I = IL 4 1 

DELP = (P (K I ) - PIKH/3.C 
P I L ) = Pi K ) 4 DELP 
Pi IL ) = P ( L ) + DELP 
S CONT IN LE 


*** THROUGH 12, EXPAND P ARRAY BY FACTOR CF 20 TC FCRM INITIAL 
EN ERGY ARRAY *** 


LMN = J 4 1 
D = 2C ,C 
K IM = 20 

KS = 1 
K IC = 3 


LOOP THROUGH STATEMENTS 26,13 ON 125 P POINTS AT A PASS 


IF(J«LT«125) GO TO 10 
KI= 125 
GU TO 11 
1C K I=J 

11 J1 = C 

K IC = K IC 4 1 

CO 12 11 = K S , K I 

DL = { F ( I 1 4 1 ) - P(I1))/D 

DO 20 KIL = 1 , K IM 

X = FL CAT I K I L ) 

XI = X-l.C 

MN = J ] 4 K IL 
2 C ESI MN ) = PUD 4- X1*DL 

12 J1 = J3 + KIM 


*** THROUGH 22, CALCULATE INITIAL TOTAL CROSS SECTION ARRAY 


NJ = MN 4 NJ 
MN = M N 4 1 
ESIMN ) = PI 1141) 
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c 

c 

c 


ENERGY LOUR 


DO 22 12 = 1 »MN 

SGS( 12 ) = O.C 

SGA S( 12 ) = C.O 
SGFSI 12 ) = G • 0 
INCJ = 1 

SR E 1= SQR T ( E S ( I 2 ) ) 

C 

C 

C NUCLIDE LOOP 

C 

C 

DO 22 IN = 1,NI 
JIG = C 
SGSP = C.C 
SGAP = C.C 
SGFP = C.C 

DEL = -2.1569179ES * RUN) * SREI 
WL = 4 .55 163 IE- 1 C / SREI 
VAL2 = $ IN ( DEL ) 

EX 1 a CMPLXt C . DE L ) 

C EXP 1 = CEXPIEX1) 

T 2 = C EXP 1 * VAL2 * WL 
ABT 2 = CA B S I T 2 ) 

SGRI = O.C 
J I = INCJ 4 NGI(IN) -1 
C 

c 

C € FACTOR LOOP 

C 

c 

CO 23 ]J = INCJ, JI 
SG APR » O.C 
SGFPR = C.C 
J IG = JIG 4 1 
NRES = NR ( IJ ) 

T 1 = C.C 
C 
C 

C RESONANCE LOOP 

C 

c 

CO 24 N2 = 1 ,NR E S 

PI = ESI 12 ) - EO(IJ,N2) 

P 2 = ( GMN ( I J , N 2 ) + GMG 1 1 J , N2) + GPF (I J ,N2) )/2 . 

DEN = CMP LX(P1,P2) 

DENSP = P 1** 2 + P 2**2 

BL AMR = WL * SREI /SQR T (ABSIEOII J , N2 ) ) ) 

SPCF = WL*BLAMR*GMN{ IJ ,N2) /DENSP 
SGAPR = GMGI I J , N2 ) * SPCF + SG A PR 
SGFPR = GM F ( I J , N2 ) * SPCF + SGFPR 
2« T 1 = ( C-MN ( I J , N 2 ) /2. / DE Nl * BLA PR+ TX 
V AL 3 = T1-T2 
VAL1 = CABS! VAL 2 ) 

SGAP = SGAP + SGAPR * GUN, JIG) 

SGFP = SGFP 4 SGFPR * GUN, JIG) 

SGSP = SGAP 4 SGFP 

23 SGRI = C( IN, JIG) * <VAL1**2 - ABT2**2) 4- SGRI 
INCJ = INCJ + 2 

SC-SII2) = COEFIIN)*! SGR I + PS1MUN) *AB T2 **24- SGS P/4 *0 ) 4- SGSII2) 
SGASI 12) = C0EFIINJ/4.C * SGAP + SGASU2) 

22 SC-FSI 12) b C0EFIINJ/4.C * SGFP + SGFSII2) 

IF! ESI 1 ) .NE.EMN ) GO TO 56 
C 

C FIRST OLTPUT POINT 

C 

NK = 1 4 NK 
ERINK ) -■ ESI 1 ) 

SGTRINK ) = SGS(l) 

SGARINK) = SGA SIX) 

SGFRINK) = SGFS(l) 

SGSRINK) = SGStl) - SGASU) - SGFSU) 


2SS 


3C5 

3C6 

3C7 


323 


33C 
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noon r> r> o o r> r> nooon n r-> o o o r>oor»oor>r>no nononooon 


*** T FRQU GH 2 6 v ELIMINATION THROUGH I N TERP C LAT I C N AND COMPARISON, 
INTERMEDIATE GENE RATI ON f I NTE RPC LAT I CN , AND COMPARISON, 

ANC OUTPUT OF SUCCESSIVE 5 CC POINT ARRAYS *** 


56 KM = C 

J 1 = MN - 2 
LEP = MN 
CIL = 3.0 

TN = 0 .C 
K A = C 

CO 26 1=1, J1 


KM = NUMBER OF SUCCESSIVE ELIMINATED POINTS 


I F( KM ) 27,26,21 


CHECK ON ENDPOINT 


21 I F( I. EG .Jl) GO TO IS 


ADVANCE I SEQUENCE 


KM = KM- 1 
GO TO 26 

IS IF( ER(NK )— ESIJ1+1) )5S,5S,26 
5S M = MN 
C-0 TO * 2 
2 6 L = 1 + 2 
DE = 2.C 
KM = C 
KJ = 2 
2 S ML = 1 + 1 
NL = L - 1 
XD = C .0 
DO 20 K=ML,NL 


INTERPOLATION TYPE, 61 SEMI LOG — 60 LINEAR 


I PC INTYP) 6 0 , 6 1 , 6C 

6 C SC = ( <SGS(L ) — SG S < I ))/(ES(L)— ES(I) ) )*<ES(K)-ESIl) ) +■ S GS I I ) 
GO TO 6 2 

63 RATIO = { E S ( L ) - ES(K))/(ES(LI - ESIH) 

SC = SC-S(L) * < SG SCI l/SGS(L))**RATI C 


COMPARE CALCULATED AND INTERPOLATED CROSS SECTIONS FOR 
ABSOLUTE AND PERCENTAGE D I FFE RE NCE S I 31 — MUST eE RETAINED) 


AC A 

62 XC = SC - SGSIK ) 

XD = XC + XC/SGS(K)*iOG.O 
DIF = ABSISC - SGSIK)) 

IF! C IF .CT .FD } GO TO 31 
PC = D 1F/SGSIK )* 1CC.0 
I F( PC • C T «PD ) GO TO 31 


CHECK FOR INFLECTION POINTS 


27 


c 


c 

c 

c 

c 

c 


c 

c 

c 

c 

c 

c 

c 

c 


c 

c 

c 


Cl = XC/ABSC XC ) 

CIM = C I* C IL 
CIL * Cl 

IF (C IM - C . 5 ) 32,33,33 


CHECK FOR MAXIMUM, MINIMUM PCIMS 


32 I F ( SG S ( K ) - SG S C I ) ) 34,34,35 

34 IF(SGS(K) - SGS(L)) 32,3C,3C 

35 i F ( SGS (K ) - SG $ { L ) ) 3C,3C,32 
3 C CONTINUE 

KM = KM +1 
CX « XC 
L = L4l 

IF { L * C T .LEP ) GO TO 31 


RETIRN FOR SUCCESSIVE ELIMINATION 


GO TO iS 

COUNTER ON MA X , M,I N , I NF LEC TI C N 

32 ICU = 1CU 4 1 
3 1 IF (KM ) 3t t 37,2* 


CHECK MAGNITUDE OF SIGMAS TC DETERMINE IMPORTANCE OF 


C PERCENTAGE DIFFERENCE AND SET UPPER LEVEL 

C 

C 

31 IF(SGSU) 4 SGS(K) + SGS(L)-2. C) 36,38,39 

38 CP = 2C.0 
GO TO K 

3 5 CP = 1C.0 

4 C IF(PC.CT.CP) GO TO 41 
C 
C 

C RETAIN POINT, GO TC STORAGE 

C 

c 

M = K 
GO TO *2 
C 

c 

c 

C 41 5 C , CALCULATE INTERMEDIATE POINTS UNTIL UPPER LIMIT 

C CRITERIA IS MET THROUGHOUT. ENERGY INTERVAL IS DIVIDED BY 

C SUCCESSIVE POWERS OF 2 UNTIL STCPED AT 256. 

C CROSS SECTION CA LC ULA TI ON ,1 N TE R PC LATI CN AND CCMPARISCN 

C SEQUENCING ARE THE SAME AS ABCVE. 

C 

c 

c 

41 KJ = 2*KJ 

I F ( KJ .LT.2SS) GO TO 43 
WRITE! 8, 1C7) E S ( I ) ,ES(L) 

GO TO ^4 

43 ET( 1) = ESU ) 

ST( 1 ) = SGS(I» 

S A T ( 1 ) = SGAS( I ) 

SFT( 1 ) = 5GF S ( I ) 

ET(KJ 4 ] ) =E S ( L ) 

S T( K J 4 ] ) = SGS(L) 

SAT(KJ41) = SGAS(U) 

SFT(KJO) = SGF 5 ( L ) 

DE = 2.0 *DE 

CL = ( ES( L )-ES< I ) ) /DE 


457 
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o o o n r> 


DO 45 N 1 = 2 , K J 
ET(N1) = E 7 < N 1— 1 ) +DL 
IN CJ = 1 

SREN I = SQR T(ETCNl)) 

4 62 

ST(N1) = C.O 
SA 7 ( N 1 ) = C.C 
SFT(N 1) = C.C 
CO 4 5 IN = 1 , N I 
J IG = C 
SGSP = C.C 
SGAP = C.C 
SGFP = C.C 

DEL = -2.1565175E5 * RUN) *SREM 
WL = 4 . 551631 E— 10/ SRE NI 
VAL2 = SIN(DEL) 

450 

EX 1 = CMPLXl C.,-DEL) 
C EXP 1 = C E XP { E X 1 ) 

4 51 

T 2 = C EXP 1 * VAL 2 * WL 
A B7 2 = CAB S ( 12 ) 

452 

SGRI = O.C 

J I = INC J 4 NGI (IN) - 1 

DO 46 IJ = INC J , J I 

SGAPR = C.C 

SGFPR = O.C 

JIG = JIG 4 1 

NR E S = NR ( IJ ) 

T 1 = C.C 

DO 4 7 N 2 = I , NR E S 

PI = E 1 ( N 1 ) - EO ( I J , N 2 ) 

P 2 = ( GMN ( I J , N 2 ) + GMG (I J * N2 ) 4 G MF (I J , lv2 ) ) /2 .0 
DEN = CMP L X ( P 1 , P 2 ) 

DENSP = P 1** 2 4 P 2** 2 

BL AMR = WL * SRENI/SQR7(ABS(E0(I J, tZ) ) ) 

5ce 

SPCF = WL * BL AMR * GMN ( I J , N2) /DE N S P 
SGAPR = GM G ( I J , N 2 I * SPCF 4 SGAPR 
SGFPR = GM F ( I J , N 2 ) * SPCF 4 SGFPR 
T 1 = ( CMN( IJ,N2)/2./DEN)*BLAMR 4 71 
VAL 3 = 11-72 
VAL 1 = CABSt VAL3 ) 

515 


SGAP = S GAP + SGAPR * GUN, JIG) 

SGFP = SGFP + SGFPR * GUN, JIG) 

SC-SP = SGAP 4 SGFP 

46 SGRI = G( IN, JIG) * <VAL1**2 - A872**2> + SGRI 
IN CJ = IN C J + 2 

ST(N1) = COE F ( IN)*(SGRI + PSTM(IN) *A 8 72**2+ SGSP/4 .0 ) +ST(M) 

SAT(Nl) - COEF(IN)/4.0 * SGAP + SA7(M) 

45 SF 7 ( N I ) = C0EF(IN)/4.C * SGFP + SF7(M) 

CO 48 N 1 = 2 , K J 

LL = N ] 4 ] 

IL = N ] - 1 

I F( IN7YP ) 64, 65,64 

64 SK = ( (STILL )- S7 { IL ) ) / (E T ( LL)-E 7 < I L) > ) * <ET { M ) -E7 ( I L > ) +ST(Itl 
GO 70 6 6 

65 RA7IG = (E7W.L) - E 7 ( N 1 ) ) / (E7 ( L L ) - E7ULI) 

SK = S 7 (L L ) * <S7(IU/ST(LL) )**RATIC 556 

66 PK = ( ABS( SK-S7(N1> ) ) / S7 ( N 1 ) *1 C C. 0 
I F{ PK .G7.DP) GO 70 41 

46 CONT IN l E 

4 4 KC = KC 4 KJ- 2 
KJ = KJ 4 1 
KQ = NK 4 KJ 
I F ( KQ .C7.5C1) GO 70 45 


0U7PU7 STORAGE ON I N7ERMEDI A 7E PCI N7S 


5C DO 51 N 1 = 2 , K J 
NK = NK 4 1 
ERINK ) = E T ( N 1 ) 

SGSR(NK) = S7CN1)— SAT(Nl) — SF7IM) 
IF( SGSRINK ) ) 67,67,68 
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67 SGSR(NK) = 0.C1 

SGTR(NK) = SAT(Nl) + SFT(Nl) + C. 01 
GO TO 6S 

66 SGTR(NK) = ST(N1) 

6S SGAR(NK) = SAT(Nl) 

EC I F = ER { NK ) - ER!NK-1> 

FLGER - ALQG ( ER ( NK ) /E R ( NK- 1 ) ) 6C 1 

ACRAKNK) = ( $GAR(NK)-SGAR (NK-1 > ) /EQI F*(EDIF-ER (NK-1) *FLGER) 

1 + SGAR(NK-1 )*FLGER + ADRAI (NK— 1) 

51 SGFR(NK) = SFT(Nl) 

KM = 1 

GO TO 26 
C 

c 

C CO ENTER ON ELIMINATED POINTS 

C 

C 

36 KA = KA + KM 
C 
C 

C CO ENTER ON PERCENT DEVIATION 

C 

C 

TN - TN + DX 
M = L- 3 
C 
C 

C IF OUTPUT ARRAY IS FULL, PRINT, RETAIN LAST PCINT, ANC EMPTY 

C 

C 

<2 I F ( NK .LE. SCI) GO TO 53 

49 WRITE (6, ICS ) 616 

WRITE! 6, ICS) (ER(KAM), SGTR(KAM), SGSP(KAN), SGAP(KAM), SGFF(KAM), 

1 ADRAI(KAM), KA M - 2,NK) 619 

NKST = NK ST + NK - 1 
I F { PC F .EQ * C . C ) GO TO 52 

WRITE! 6, 132)(ER(KAM),SGTR(KAM) , SG SR ( K A N ) , SG AR ( K A N) ,SGFR(KAM), 

1 KAM = 2, NK ) 622 

52 IF! E$(M ).EQ.EM>) GO TO 77 
ER! 1) = ER (NK ) 

S GAR ( I ) = SGAR (NK ) 

ADRAI! ] ) = ADRA I ( NK ) 

CO S4 NK S = 2 , NK 
ER ( NK S ) = C.C 
SGTR(NKS) = C.C 
SGSR(NKS) = C.C 
SGAR! N KS ) = C.C 
54 SGFR(NKS) = C.C 
NK = 1 
C 

c 

c RET CRN TO INTERMEDIATE PCINT SEQUENCING 

C 

C 

IFCKQ .C-T.5C0 ) GO TO 50 
C 
C 

C REGELAR STORAGE OF RETAINED FCINTS 

C 

C 

52 NK = NK + 1 

ER(NK ) = E S ( M ) 

C 

c 

C 

C CORRECTION FOR NEGATIVE SCATTERING CRCSS SECTION ARISING 

C FROM EXCESSIVE INTERFERENCE DEE TC INADEQUACY OF G F ACT GR= 1/2 

C APPROXIMATION AT S hAVE DIP BELCF SMALL RESONANCE SFACCWEC 

C BY LARGE RESONANCE RARE SITEATICN OCCURRING ONLY WHEN SPIN 

C STATES ARE UNKNOWN 

C 

C 
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onnon ononn nnooo r> n o r> o n r>r> ri r> n r> r> 


S GSR < NK ) = SGS(M) - SGAS(M) -SGFS(M) 
IFtSGSR(NK)) 7C,7C,n 
7 C SGSR( NK 1=0.01 

SGTR(NK) = SGA$(M) 4 SGFS(M) 4 0. 01 
GO TO 72 

71 SC-TR(NK ) = SGS(M) 

72 SGAR(NK) = SGA S<M ) 

ECIF = ER ( NK I - ER < NK- 1 ) 

FLGER = AL0G(ER<NK)/ER<NK-1)) 


ACCUMULATIVE OILUTE RESONANCE ABSCRFTICN INTEGRAL 


ADRAUNK) = ( SGAR(NK)— SGAR(NK— 1) )/E0IF*(EDl F— £R (NK— 1 ) *FLGER ) 
1 4 SGAR(NK-1)*FLGER4 ADRAI(NK-l) 

SGFR(NK) = SGFS(M) 

2 6 CONTINUE 

COUNTERS ON KA,TN 

KOI = KOI 4 KA 
CKO 2 = CKO 2 4 TN 


LOOP ON P POINTS 


K S = K ] 4 1 

K CM = KIC * 3 2 5 
IF(J.LT.KCM) GO TO 12 
K I = K CM 
GO TO 11 

12 IF(K I .EQ .J ) GO TO 14 
K I = J 
C-0 TO 3 1 


LOOP ON RESONANCES 


14 ND = MC 4 1 

I CM = JC * 8 2 
IFtNDES.LT ■ ICM J GO TO 18 
MC = ICM 
GO TO ] 6 

16 IF(MD.EQ.NDES) GO TO 76 
MO = NCES 
GO TO 36 


FINAL OUTPUT 


76 WR ITE ( 6, 1C9 ) 721 

WRITE(6,1C5) (ER(KAM), SGTR(KAP), SGSR(KAM, SG/P(KAM), SGFP(KAM), 

1 ACRAI(KAM), KAM = 2,NK) 722 

NKST = NK ST 4 NK - 1 
IFIPCF.EQ .C .0) GO TO 77 

WRITE! 6 , 1 32 HER (KAMI , SGTR ( KAMI ,SGSR(KAP> ,SGAR ( K AM » SGFR( KAMI, 

1 KAM = 2 y NK ) 727 

7 7 TN = CKO 2 
KA = KOI 
NJ = NJ 4 1 


SUMMARY OF CALCULATION 


WR ITEC 1C6) N f NKST,FD f PD,NJ,KA f KC,ICU 750 

ON = FLOAT(KA) 

TN = TN/DN 

WR ITE( l« 1 C E ) TN 751 
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ICC FORMAT ( 1H 1 ) 

1C j F0RMAT(2I1C,4F1C.5,3I5,F5.C> 

10 2 FORMAT (2F 1 C . 5 , E 1 1 . 5 , I 9 ,F 1 C. 5 ) 

103 F0RMAT(F1G.5, I 1C) 

104 F0RMAT(F1C .5 , I 1C, 3F10.5) 

105 FURMAT ( 1HC6F2C.9) 

106 FORMAT ( 1HC 2CX 2HFQR 15, 2X11H RE SO NANCES, 16, 2 X54HE N E PGY POINTS ARE REQ 
1U IR ED FOR DEVIATIONS OF LESS TH AN/ /20 XF5 . 1 ,2 X9 HB ARNS A NCF5 * 1, 2X 3 1H 
2PERCENT AT INTERPOLATED POI NTS. //2 0 >14HCF THE I MT I AL 15 ,2X7FP0 INTS 

3, 

415, 2X16FWERE DISCARDED AND I 5 ,2 X32HADDI T I CNAL POINTS WERE INSERTED. 
5//2CX If, 2X46HPGINTS WERE RETAINED AS PAXfMIN CR INFLECTION.////) 

107 FORMAT {1HC2CX14FREGI0N BE TWEENF1 0. 5 ,2 X6HE V ANDF10.5 , 2X48HEV CANNOT 
1 BE FITTED WITH 256 INTERMEDIATE POINTS//) 

10e FORMAT! 1H C // 20 X 45HA VERAGE PERCENT DEVIATION GF ClSCARDEC PC INTS// 4 
15XF10.5 ) 

109 FORMAT < 1HC6X6 HE NERGY14X11H SIGMA TCTAL9X13HSIGMA SC ATT E P7X 13hS I CM A 
1CAPTURE7X13HSIGMA F I S SI 0N6 XI 4H ACC. I (DILUTE)//) 

UC FORM AT ! 1HL 1 IX E FMA S S N0.3X8HPCT. AB.5X11HNUC. R AC IUS3X7 HG V ALU E5X 1 1 
1FRES. ENERGY2X7HGAMMA N5 X 7HGAMNA G5X7HGAMMA F//) 

111 FORM AT ( 1HC 1C XF 6 • 1 ,6XF 6. 4, 6X1 PE 12.6 ,3 X0FF6 .4 ,3X4F12 .5/ ) 

112 FORM A T ( 1 H$ , IP 5E 12.5) 

113 FORMA T ( II, IX, 14A5) 

114 FORMAT (2X, 14A5 ) 

115 FORMAT ( IHL11X34HLINEAR INTERPOLATION WILL BE USED.) 

116 FORMAT! 1HL11X43FSEMILOGARITHMIC I N TE R FC LAT I CN WILL EE USED.) 

RETURN 

ENC 
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N/TURAL TUNGSTEN 


ALL PARAMETERS, WITH THE EXCEPTIONS NC TED BE LCW , ARE TAKEN FRCM 
REFERENCE A-l, WHICH OOE S NOT LIST THE LOW ENERGY DATA. THE REMAINING 
PARAMETERS WERE CHOSEN ON THE BASES CF CONSISTENCY WITH THE MAIN 
SOURCE, MINIMUM EXPERIMENTAL UNCERTAINTY, AND BEST PARAMETERS FOR 
MATCHING MEASURED VALLES OF THERMAL CROSS SECTIONS AND DILUTE 
RESONANCE INTEGRALS. 

W-1EC PARAMETERS FROM REFERENCE A-2 . 

W-1E21-32.999EV >AND W- 1 £4 (-9 9. 97 8 EV) HAVE ARBITRARY ENERGIES C ET A IN ED 
USING THE AVERAGED REDUCED NEUTRON WIDTH AND ASSUMED CAPTURE WICTH OF 
REFERENCE A-l TO MAKE LP THE DIFFERENCE BETWEEN THE THERMAL CROSS 
SECTIONS OF REFERENCE A-3 AND THE VALUES CALCULATED FRCM THE KNOWN 
PARAMETERS. 

W- 1 £2( 4.16 AND 21. C9 EVI AND W-ie3<7.67 EV! FRCM REFERENCE A- 4. 
W-1E3(27.06 EV) FROM REFERENCE A-5. 

W-1E3( 46.2, <7.£, AND 65.3 EV) FRCM REFERENCE A-6. 

W-l€3( G FACTOR * 1/2 RESONANCES) FRCM RECLCED NEUTRON WICTHS, CAPTURE 
WIDTHS, AND ENERGIES GIVEN IN REFERENCE A-l. 

W-1E6( 18.82 EV) PARAMETERS FRCM REFERENCE A-7. 


THERMAL (C.C252 EV) ABSORPTION CROSS SECTIONS (BARNS) 


W- 1 EC 

CAL C LLA T ED 
3.5 

REF. A-3 

W-162 

20 .76 

2 C . 7 ( C. 5) 

W-162 

8 .69 

1 C . C ( C. 3 ) 

W-l 64 

1.71 

1. 7(0.1) 

W-1E6 

35 . a 

37.6(1.2) 

NAT . W 

— 

18.3(0. 5) 


INFINITELY OILLTE RESONANCE INTEGRALS (BARNS) ABOVE 0.5 EV 

CALCLLA TED 
W- 1 EC 204 
W-162 623 

W-163 355 

W-164 15. E 

W-1E6 457 


CALCULATED VALLES HAVE AN UPPER LIMIT CF 2200 EV 
REFERENCES 

A-l BLOCK, R.C., HOCKENBURY ,R. W. , RL S SE L L , J. E. , THE PARAMETERS CF THE 
NEUTRON RESONANCES IN W- 182 , W-l 83 , W-l 84 , AND W-186, CRNL-3924, 

31-35, (MAY 1566) (SEE ALSO RPI-328-56). 

A-2 JUNG,H., BLOCK, R.C., SLA LGH TER ,G .G . , PARAMETERS CF NEUTRON 
RESONANCES IN W-1B0, ORNL-3924, P. 30, (MAY 1966). 

A-3 FR IESENHAHN, S.J., HADDAD, E. , F R CHNE P »F . H. , LOPEZ, W.M., THE 

NEUTRON CAPTURE CROSS SECTICN CF THE TUNGSTEN ISCTCPES FROM 0.01 
TO 10 ELECTRON VOLTS, NUCLEAR SCIENCE AND ENGINEERING 26, 

4 6 7—49 9 ( 1966 ) . 

A-4 HARVEY, J. A., THE MEASUREMENT OF TOTAL NEUTRON CRCSS SECTIONS IN 
THE RESONANCE ENERGY REGION AND THE DETERMINATION CF RACIATIGN 
WIDTHS OF RESONANCES, ANS TOPICAL MEETING CN REACTOR PHYSICS IN 
THE RESONANCE AND THERMAL REGICN.P 10, SAN DIEGC (FEBRUARY 1966). 

A-5 P A Y A , C ■ , PEARCE, K.D., HAR VE Y , J. A. , SLA UGHTE R ,G .G. , PARAMETERS OF 
LOW ENERGY RESONANCES IN TUNGSTEN, CRNL-3582 ,58-60 (JUNE 1964). 

A-6 RUSSELL, J .E HOCKENBURY ,R. W. , BLCCM.R.C., NEUTPCN CAPTURE 
MEASUREMENTS ON THE ISOTOPES OF TUNGSTEN, WASH-1046, P. 104 
(JANUARY 1964). 

A-7 PEARCE, C.R., SHOOK, D.F., DETERMINATION CF TUNGSTEN RESONANCE 

ABSORPTION INTEGRALS BY ACTIVATION, NASA TN D- , (1967). 


REF. A-7 

15(2) 

440(20) 



LINEAR INTERPOLATION hlLL BE LSED 


MASS NJ. 

PC T . AB. 

NUC. RADIUS 

G VALUE 

RES. ENERGY 

GAMMA N 

GAMMA G 

GAMMA F 

i ec.c 

C.0C13 

E. 4 7C50GE-13 

1.0000 

15.90000 

0 • C 1 2 f C 

C.C60C0 

-0. 

lec.c 

C.C013 

8. 47050CE-13 

1.0000 

49.30000 

0 .C062C 

0.C60C0 

-O'. 

1 tc .c 

C.CC13 

e. 47050CE-13 

1.0000 

62.70000 

0 .CO 1*0 

O.C6CCO 

-0. 

lec.c 

G.OC13 

8. 47G50GE— 13 

1.0000 

75.20000 

0.035CC 

O.C6CCO 

-0. 

itc.c 

C.0C13 

8.470500E-13 

1.0000 

87.40000 

0.CC72C 

0.C6CC0 

-0. 

MASS NO. 

PC T • AB. 

NUC. RADIUS 

G VALUE 

RES. ENERGY 

GAMMA N 

GAMMA G 

GAMMA F 

1 £2 . C 

C.2631 

8. 5CC590E-13 

1.0000 

-32.99900 

0. 1C514 

C.C7CC0 

-0. 

1 £2 . C 

C . 263 1 

8. 5 COS 9 0E-13 

1.0000 

4.16000 

C. 00148 

0.C5400 

-0. 

1 £ 2 • C 

C.2631 

8. 5CC590E-13 

1.0000 

21.09000 

O.C423C 

C.C62C0 

-0. 

1 £2 . C 

C . 2631 

8. 5C0590E-13 

1.0000 

98.00000 

0.000C1 

C.C7CC0 

-0. 

1 £ 2 . C 

C.2631 

B. 5 C0590E-13 

1.0000 

114.70000 

C.290C0 

0 .C65C0 

-0. 

1 £2 . C 

C.2631 

8. 5CC59CE-13 

1.0000 

130 .50000 

O.CC050 

0.C7OCO 

-0. 

1 £2 . C 

C.2631 

8. 5 CC590E-13 

1.0000 

214.00000 

0 .00 2 £ C 

C.C70C0 

-0. 

1 £ 2 • C 

C.2631 

8. 5 CC59QE-13 

1.0000 

250.00000 

1.100C0 

C.C64C0 

-0. 

1 £ 2 . C 

C.2631 

8. 5 CC590E— 13 

1.0000 

283 .00000 

0.CC11C 

C .C70C0 

-0. 

1 £2 . C 

C.2631 

8. 5CC590E— 13 

1.0000 

303.00000 

0.00015 

0.C7CC0 

-0. 

162. C 

C.2631 

8. 5C0590E-13 

1.0000 

343.00000 

o .cce*c 

O.C7CCO 

-0. 

1 £ 2 • C 

C.2631 

8. 5 CC590E-13 

1.0000 

378 .00000 

0.130CC 

0.C65C0 

-0. 

1 £ 2 . C 

C.2631 

8. 5 C0590E — 13 

1.0000 

410 .00000 

0 .C002C 

0.C70C0 

-0. 

1 £2 • C 

C.2631 

8. 5 CG590E-13 

1.0000 

430.00000 

C.280C0 

0.C55C0 

-o. 

1 £2 . C 

C.2631 

8. 5 CC590E-13 

1.0000 

486.00000 

0.500CC 

0 .C49C0 

-0. 

1 £2.C 

C.2631 

8. 5 CC590E— 13 

1.0000 

580.00000 

C.3C0C0 

O.C57CO 

-0. 

1 £ 2 . C 

C.2631 

8. 5CC59CE-13 

1.0000 

616.00000 

O.CC 14C 

C • C 7CC 0 

-0. 

1 £ 2 • C 

C. 2631 

8. 5 CC590E-13 

1.0000 

658.00000 

0. 160CO 

0 .C5900 

-0. 

1 £ 2 . C 

C.2631 

8. 5C0590E-13 

1.0000 

673.00000 

0.C0060 

C.C70CO 

-0. 

1 £2 . C 

C.2631 

8. 5CC59CE-13 

1.0000 

762.00000 

0.069CC 

0 .C74C0 

-0. 

1 £2 . C 

C.2631 

8. 5 CCS 90E-13 

1..0000 

785.00000 

C.C14CC 

C.C70C0 

-0. 

1 £2 . 0 

C.2631 

8. 5CC59CE-13 

1.0000 

866.00000 

0.C21C0 

C.C7CC0 

-o. 

1 £ 2 • C 

C.2631 

8. 5 C059 0E-1 3 

1.0000 

922.00000 

C . 4000 C 

0.C78C0 

-o. 

1 £2 . C 

C.2631 

8. 5CC590E-13 

1.0000 

951.00000 

2.200CC 

0 .C77C0 

-0. 

1 £ 2 . C 

C.2631 

8. 5C0590E-13 

1.0000 

1010.00000 

0.45G0C 

0.C7CC0 

-0. 

1 £2.C 

C.2631 

8. 5 CC590E-13 

1.0000 

1100.00000 

1.600CO 

C.C65C0 

-0. 

1 £2 * C 

0.2631 

8. 5 CC590E-13 

1.0000 

1170.00000 

0.480CC 

O.C58CO 

-o. 

MASS NO. 

PC T • AB. 

NLC . RADIUS 

G VALUE 

RES. ENERGY 

GAMMA N 

GAMMA G 

GAMMA F 

1 6 2 . C 

C. 142 6 

£• 5 161 2 OE-13 

0.2500 

47.80000 

C.l 15CC 

C.C76C0 

-o. 

1 £3 . C 

C. 1426 

6. 51612 CE— 1 3 

0.2500 

144.60000 

0.C95C0 

C.CSOCO 

-o. 

1 63 . C 

C.1428 

£. 516120E-13 

0.2500 

154.90000 

C .4 IOC C 

o.cecco 

-0. 

163. C 

C. 1426 

8. 51612 OE-13 

0.7500 

7.67000 

O.C0174 

C.C79C0 

-0. 

163. C 

C.1428 

6.51612CE-13 

0.7500 

27.08000 

0 .04 22 C 

C.C75C0 

-0. 
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183.0 

0. 142 8 

8. 51612 QE-1 3 

0.7500 

46.20000 

0.154CC 

0 .C68C0 

-0. 

1 83. C 

C . 14 2 8 

8. 5 1 612 OE-13 

0.7500 

65.30000 

C.CC18C 

O.C75CO 

-0. 

18 3.0 

C . 1428 

8. 51612 CE— 13 

0. 7500 

10 1.300 00 

o.cseco 

0 .CS8C0 

-0. 

183. C 

C. 1428 

8. 51612 CE-13 

0.7500 

157.30000 

0 .C67CC 

C.C80C0 

-0. 

183. C 

C.1428 

8. 51612 OE-13 

0. 7500 

174.30000 

0.C53C0 

0 .csoco 

-0. 

183. C 

C.1428 

8. 51612CE-13 

0. 7500 

192.60000 

0 .C25CG 

O.C7CCO 

-o. 

183.0 

C. 1428 

6. 51612 OE-13 

0.7500 

259.00000 

0.C50C0 

O.ICOCO 

-0. 

183. C 

C.1428 

8. 5 1612 CE-13 

0.7500 

280 .00000 

0.21CCC 

0 .C73CO 

-0. 

183. C 

C. 1428 

8. 5 161 2 OE-13 

0.7500 

297.00000 

O.O29C0 

C.C99C0 

-0. 

183. C 

C.1428 

8. 516120E-13 

0.7500 

323.00000 

0 . C 8 2C C 

C.C96C0 

-0. 

183. C 

C.1428 

8. 516120E-13 

0.7500 

349.00000 

0.140C0 

O.C95CO 

-0. 

183. C 

C.1428 

6. 516120E-13 

0. 7500 

361.00000 

0.C35CC 

0 .C7CC0 

-0. 

183. C 

C.1428 

6. 51612 OE— 13 

0.7500 

379.00000 

0.C740C 

0.C78C0 

-0. 

MASS NO. 

PC T . AB. 

NUC. RAC I US 

G VALUE 

RES. ENERGY 

GAMMA N 

GAMMA G 

GAMMA F 

183. C 

C.1428 

8. 5 1 61 2 CE-13 

0.5000 

40 .70000 

0 .0C6 IS 

C.C80C0 

-0. 

183. C 

C.1428 

6. 516120E— 13 

0.5000 

104.10000 

0.00619 

O.C0OOO 

-0. 

1 83. C 

C.1428 

8. 51 61 2 OE-13 

0.5000 

138.30000 

0.CC564 

c .ceoco 

-0. 

183. C 

C.1428 

E. 51612 OE-13 

0.5000 

203.80000 

0.00 3*3 

c .cecco 

-0. 

183. C 

0.1428 

8. 51612 CE-13 

0.5000 

220.80000 

0.00525 

c .cecco 

-0. 

183. C 

C. 1428 

8. 516120E-13 

0.5000 

228 .00000 

O.OClf 1 

c. ceoco 

-0. 

183. C 

C.1428 

8. 516120E-13 

0.5000 

236.00000 

0.01643 

C.C80C0 

-0. 

183. C 

C. 1426 

8.51612 0 E— 1 3 

0.5000 

241.00000 

0.05278 

o .ceoco 

-0. 

1 83. C 

C. 1428 

8. 51612 OE— 13 

0.5000 

244.00000 

O.CC553 

c .cecco 

-0. 

183. C 

C.1428 

8. 51612 OE-13 

0.5000 

289 .00000 

0 .CO 7 * 8 

0 .C8CC0 

-0. 

183. C 

C.1428 

8. 516120E-13 

0.5000 

338 .00000 

0 .C 20 20 

c. ceoco 

-0. 

183. C 

C.1428 

8. 516120E-13 

0.5000 

354.00000 

0 .CC53C 

c .cecco 

-0. 

183.0 

C.1428 

6. 51612 CE-13 

0.5000 

392 .00000 

C .C475C 

o. cecco 

-0. 

MASS NO. 

PC T . AB . 

NIC. RADI IS 

G VALUE 

RES. ENERGY 

GAMMA N 

GAMMA G 

GAMMA F 

18 4.0 

0.3064 

8. 53166CE-13 

1.0000 

-99.97500 

0.24650 

0.C70C0 

-0. 

184. C 

C. 3C64 

8. 53166CE-13 

1.0000 

102. 10000 

C.CC35C 

0.C7CC0 

-0. 

184. C 

C.3C64 

8. 53166CE-13 

1.0000 

184.70000 

1.20000 

O.C78CO 

- 0 . 

184 .C 

C.3C64 

8. 5316 60E— 13 

1.0000 

244.00000 

0.CC24C 

C .C7CC0 

- 0 . 

184. C 

C • 3064 

8. 531660E-13 

1.0000 

311.00000 

0.075CC 

0 .C51C0 

-o. 

184. C 

C.3C64 

8. 53166CE-13 

1.0000 

424.00000 

0.042CO 

0.C69C0 

-0. 

184. C 

C.3C64 

8. 53166 CE-13 

1.0000 

595.00000 

0 .C01C0 

C .07CC0 

-0. 

184. C 

C.3C64 

8. 53166CE-13 

1.0000 

604.00000 

0 .68OC0 

0 .C78CO 

-0. 

184.0 

C • 30 64 

8. 5316 6CE— 13 

1.0000 

705.00000 

0.CC660 

C.C7CC0 

-0. 

184.C 

C. 3C64 

8. 53166CE-13 

1.0000 

787.00000 

O.C260C 

C.C70CO 

-0. 

184. C 

C. 3C64 

8. 5 31 66CE-13 

1.0000 

802.00000 

1.6C0C0 

C.C7SC0 

- 0 . 

184.0 

C • 3064 

8. 5316 60E— 13 

1.0000 

961.00000 

1.60000 

0.11CC0 

-0. 

1 8* . C 

C • 3064 

8. 5 3166 OE— 13 

1.0000 

1000.00000 

0.140CC 

C.C55C0 

- 0 . 

184. C 

C . 3064 

8. 53166CE— 13 

1.0000 

1090.00000 

3 .40000 

0.CS5C0 

-0. 
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1E4.C 

C.3C64 

8. 53166GE-13 

1.0000 

1140.00000 

0.340CC 

0.C6CC0 

-0. 

1E4.C 

C . 3064 

8. 53166CE-13 

1.0000 

1270.00000 

1.2G0CC 

0.C63C0 

-0. 

1E4 .C 

C.3C64 

8. 5 31 66CE-13 

1.0000 

1340 .00000 

0.CC21C 

C.C70C0 

-0. 

iE4.o 

C.3C64 

8. 531660E— 13 

1.0000 

1410.00000 

2.7C0CC 

0.C76C0 

-0. 

1E4.C 

C.3C64 

8. 531660E-13 

1.0000 

1430.00000 

C.250CC 

C.C75C0 

-0. 

1E4.C 

C • 3064 

8. 53166CE-13 

1.0000 

1520.00000 

1.300CC 

0.C58C0 

-0. 

IE* .0 

C.3C64 

8. 53166CE-13 

1.0000 

1560.00000 

O.C66CC 

0.C7CC0 

-0. 

164.0 

C • 3064 

8. 531660E-13 

1.0000 

1660 .00000 

0.25000 

o.ceoco 

-0. 

164.0 

C.3C64 

8.531660E-13 

1.0000 

1800.00000 

1.1C0GC 

C.C65C0 

-0. 

164.C 

C.3064 

8. 531660E-13 

1.0000 

1880.00000 

0.C3CC0 

0.C7CC0 

-0. 

1 £*.C 

C.3C64 

8. 531 66CE-13 

1.0000 

1930.00000 

C.250C0 

0 .C75CO 

-o. 

MASS NO. 

PC T . AB. 

NIC. RADIUS 

G VALUE 

RES. ENERGY 

C-AMHA N 

GAMMA G 

GAMMA F 

1 66 . C 

C.2664 

6. 56242 OE-13 

1.0000 

18.83000 

0.315CG 

C.C41C0 

-0. 

166. C 

^ • 2864 

8. 56242 OE-13 

1.0000 

111 .30000 

0.00004 

0 .C 7000 

-0. 

1 U .C 

C.2E6* 

6. 56242 OE— 13 

1.0000 

171.50000 

0.C27CC 

C.C75C0 

-0. 

166 . C 

C.2E64 

8. 56242 OE-13 

1.0000 

197.60000 

c.cccec 

C.C7CC0 

-0. 

166. C 

C.2864 

8. 56242 OE— 13 

1.0000 

218 .00000 

C.530CC 

0.C62C0 

-0. 

166. C 

C .2864 

8. 56242 C E-13 

1.0000 

245.00000 

0 .00020 

C.C7CC0 

-0. 

166. C 

C.2664 

8. 56242 OE-13 

1.0000 

288 .00000 

0.C26CC 

0.C7CC0 

-0. 

1 66. C 

C.2664 

8. 56242 OE— 13 

1.0000 

407.00000 

O.C75CC 

C.C62C0 

- 0 . 

166. C 

G . 2864 

8. 56242 OE— 13 

1.0000 

458.00000 

o.cooeo 

C.C70C0 

-0. 

166. C 

C.2664 

8. 56242 OE — 13 

1.0000 

512.00000 

0.C65CC 

0.C45C0 

-0. 

166. C 

C.2664 

8. 56242 OE-13 

1.0000 

543.00000 

C.5C0C0 

C .C65C0 

-0. 

1 66. C 

C.2864 

8. 562420E— 13 

1.0000 

666.00000 

0 • 7 500 C 

o.cesco 

-0. 

166. C 

C. 2664 

8. 56242 OE-13 

1.0000 

732.00000 

2.1CCCC 

C.C52C0 

-0. 

166. C 

C.2664 

8. 56242 OE— 13 

1.0000 

774.00000 

C.0060C 

0.C7CC0 

-0. 

166. C 

C.2864 

8. 56242 OE— 13 

1.0000 

835.00000 

0.C2CCC 

C.C7CC0 

-0. 

1 66. C 

C.2864 

8. 56242 OE-13 

1.0000 

858 .00000 

0.CC2CC 

C.C7CC0 

-0. 

166. C 

C.2664 

6. 56242 OE-13 

1.0000 

968 .00000 

1. 1C0CC 

C.C72CO 

-c. 

166. C 

C.2664 

8. 562 42 OE— 13 

1,0000 

1080.00000 

0.65CC0 

O.CfilCO 

-0. 

166. C 

C.2664 

8. 56242 OE-13 

1.0000 

1130.00000 

C.45CCC 

0 .C62C0 

-0. 

166. C 

C.2664 

8. 56242 OE-13 

1.0000 

1190.00000 

0.770CC 

0.C55C0 

-0. 

1 66. C 

C.2664 

8. 56242 OE-13 

1.0000 

1420.00000 

C.250CO 

0.C50C0 

-0. 

1 66. C 

C.2864 

8. 56242 OE-13 

1.0000 

1510.00000 

1.20000 

0 .C68C0 

-0. 

1 66. C 

C.2664 

8. 56242 OE-13 

1.0000 

1550.00000 

0.C046C 

0.C70C0 

-0. 

166. C 

0.2864 

8. 56242 OE— 13 

1.0000 

1800.00000 

C.100C0 

C.C65C0 

-0. 

1 66. C 

C.2864 

8. 56242 OE-13 

1.0000 

1940.00000 

0 .55CCC 

0.C6CCC 

-0. 

166 .C 

C.2664 

8. 562 42 OE— 13 

1.0000 

2040.00000 

0.4C0CC 

0.C67C0 

-0. 

166. C 

C.2664 

e. 56242 CE— 1 3 

1.0000 

2120.00000 

O.llOCO 

C.C67C0 

-0. 

*CY 

SIGMA TOTAL 

SIGMA SCATTER S I GM A 

CAPTURE 

SIGMA FISSION ACC. I <DI LUTE1 

C2 5255556 

21. 

5686655C4 

4.110808134 

1 7 .457 857 3 7C 

-c . 

-0. 000000000 

0^2255747 

17. 

644255 S3 0 

4. 109911680 

13.534348245 

-C. 

8.050632238 
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O.CfS2S54S9 

15.567874193 

4.109017253 

11.458856540 

-C . 

12.291931272 

0 • C “16299243 

14.23493 862 2 

4. 108123541 

10.126815081 

-c. 

15. 019442558 

C .C 9229 £994 

13.267731886 

4.107230544 

9.180 501242 

-c. 

16. 964428186 

0.1 1C296745 

12. 57C75 6462 

4.106336263 

8.464420159 

-c . 

18.442830801 

0.1^7296497 

12.004122496 

4.105446517 

7.890675578 

-c. 

19.616559505 

0. 1*4296241 

11. 542C1 63 87 

4.104555130 

7.437461257 

-c. 

20. 578337193 

0.161297992 

11.1 56C49490 

4.103664279 

7.052385211 

-c. 

21.385611057 

C. 176297743 

1C. 827566385 

4.102773666 

6.724752715 

-c. 

22.076130629 

0.195297495 

10. 543737650 

4. 101882994 

6.441854656 

-c. 

22. 675863504 

0.212297238 

10.295413613 

4. 100992799 

6.194420815 

-0. 

23.203340054 

0.2 29296 99C 

10. 075668607 

4.100102782 

5.975765624 

-c. 

23.672186851 

0.246296741 

9. 88CC35400 

4.099212706 

5.7808226S4 

-0. 

24. 092679501 

0.263296492 

9. 7C4C12 871 

4.098322690 

5.605650161 

-c. 

24.472733021 

C.26C296236 

9.544743180 

4.097432613 

5.44731056? 

-c. 

24. 818557262 

0.297295995 

9.395791241 

4.096542478 

5.303248763 

-c. 

25.135101795 

C. 6*1267593 

7. 53 C240344 

4.067599654 

3 .46264 666C 

-c. 

29. 857317209 

1.792272965 

7. CC 5766161 

4.006211638 

3.003554452 

-c. 

32.326045036 

2.8 756 36CC 6 

8. 1 405 0C416 

3.847042620 

4.293537796 

-c. 

34. 025119781 

2.266772CC2 

10.235445510 

3.695003390 

6 .540446 5 2C 

-c. 

34.712892056 

2.527526C16 

14.350287318 

3.505337715 

10.8 44945603 

-c . 

35. 378332615 

2.6579C5996 

15.220681429 

3.361116529 

15.8595645CC 

- c , 

35. 862383842 

2. 723C950CC 

23.451154709 

3.273392439 

20.177762270 

-c . 

36. 180564404 

2. 7682 64CC4 

3C.055C46469 

3.178499460 

26 .9 2054 7C05 

-c. 

36.589158058 

3.820876506 

34.945631981 

3. 132136345 

31.813495636 

-c. 

36. 840725899 

2.6 f 2473CC8 

41.411273556 

3.091688633 

38.319585323 

-c . 

3 7. 138567448 

3. 6(60 675 1C 

5C. 3CC837C40 

3.065764904 

47.235072136 

-c. 

37.498832703 

2.9 16662C12 

62. 965734463 

3.071518421 

59*914216042 

-c. 

37.946254253 

3.924959263 

71.459150314 

3-096199989 

68.362950325 

-0. 

38.212457657 

3.9 512 56514 

81.955331802 

3.144636154 

78.810695646 

-c . 

38. 516617298 

3.967553765 

55. 166656494 

3.228197098 

91.938455356 

-c. 

38. 868031502 

2.963851016 

112. 105991364 

3.364455223 

108 .74 1536 14C 

-c . 

39.279355049 

4 .OCC 146237 

134.30C1 88065 

3.581350327 

130.71883 7736 

-c . 

39. 768147945 

4.C16445456 

164.133132935 

3.924852371 

160.208280563 

-c . 

40. 35963 7260 

4.022742679 

205.482116699 

4.473558426 

201.008556273 

-c . 

41. 091009617 

4.C4904C02C 

2 64. 96673 5655 

5.368614197 

259.598125456 

-c. 

42.019774437 

4.C69645C21 

387.312646773 

7.466571808 

379.846076565 

-c. 

43.65 8286571 

4.C J2226022 

5C5. 205785461 

9.702808380 

495 .506977C 6 1 

-c. 

44. 999054909 

4.094811022 

683.232276916 

13.341011047 

669.891265865 

-c. 

46. 778362751 

4.1C2122C23 

855. 62C7122 80 

17.078094482 

838.542617758 

-c . 

48.309873104 

4.1C7294023 

965. 13 856C322 

19.537063599 

945.601516724 

-c. 

49.214548111 

4.111455C23 

1094.768676758 

22.519088745 

10 72.249 58 8C 13 

-c. 

50.236818790 

4. 1 15616C24 

1248.874603271 

26.152420044 

1222.7 2216 32 28 

-c. 

51.398771286 

4.1 19777024 

1432. 64347 83 94- 

30.595443726 

1402.048034666 

-c. 

52.72 72 95876 
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4.1 23928024 

1651.987335205 

36.038360596 

1615.948974605 

-c. 

54.251762390 

4. 12SC^^C2A 

1913.112152100 

42.697601318 

1870 *4 14550761 

-C. 

56. 013356209 

4.1336*7025 

2335.235166260 

53.845367432 

2281.393798628 

-C. 

58. 804238796 

4. 1278C£025 

2 70S. C3 6495023 

64.105133057 

2644.93136556? 

-c. 

61.287439823 

4.1*6120025 

3554.989532471 

88.895812988 

3466.093715462 

-c. 

67. 426372528 

4 .l‘C251C26 

3555. 5785 82 764 

101.929351807 

3857.64923095? 

-C . 

71. 102688709 

4.1 ‘3065026 

41 63. 887207 C31 

109.858734131 

40 74 .0 2 8 47 290 C 

-C. 

73.755782127 

4.155829026 

4348.421651895 

116.502502441 

4231.919189452 

-c. 

76. 531302452 

4.157226C26 

4402.722595215 

119. 165344238 

4283 .5 57 250 57 7 

-C. 

77. 952976227 

A. 1566 12026 

4436.175592676 

121.309265137 

4314.870727535 

-C. 

79.388350487 

4.16C0CC026 

4447.822814941 

122.890380859 

4324.932424082 

-C. 

80. 829081535 

A, 161386567 

4437. 315C02 441 

123.881103516 

4313.433896926 

-C . 

82.268830299 

162773967 

4404.975952148 

124.270751953 

4280.70520C155 

-c. 

83.699986458 

KS547S67 

4279. 189086514 

123.294677734 

4155.894405160 

-c. 

86. 509539604 

4. 16622 1967 

4C84. 234405518 

120.210144043 

3964. 024 26 1475 

-c. 

89.210590363 

4. 172462567 

37C4. 152313232 

112.599822998 

3591. 55249 0234 

-c . 

92. 974343300 

4 . 1 76C 2 C 5 6 E 

3133. 866286133 

99.504547119 

3034.363735014 

-c. 

97.375096321 

4.1E2578568 

255C. 453133545 

85.949829 102 

2504.543304442 

-c. 

101. 045792580 

4.1 (7729966 

2232. 804 Cl 61 1 3 

76.580963135 

2156.223052575 

-c . 

1 03. 358903885 

4. 1519CC969 

1924.36330C781 

68.235794067 

1856.147506714 

-c. 

105.349685669 

4.196061969 

1662. 9e922 7295 

60.965011597 

1602.024215658 

-c. 

107. 064538956 

4.2CC222969 

1443.264759521 

54.703170776 

1388.581586745 

-c. 

108. 546264648 

4.2C4363969 

1255.1 C5466553 

49.337020874 

1209.772445675 

-c. 

109. 832181931 

4.2C854457C 

1 1C4. 580322266 

44. 742370605 

1059.83795166C 

-c. 

110. 953846931 

4.2127C557C 

574.528137207 

40.801406860 

933.726730247 

-c. 

111. 93 7504768 

4.221C2797C 

771.211875516 

34.479545593 

736 .73233C 22 2 

-C . 

113.585330009 

4.229349971 

622. 543305 969 

29.715682983 

593.227622566 

-c. 

114.895122528 

4.2*1822572 

467.632614136 

24.538002014 

443.09461212 2 

-c . 

116. 421575546 

4.2 ‘6476972 

335. 9CC8522C3 

19.931400299 

315 .96945 19C 4 

-c . 

117.907953262 

4.270959973 

27C.35633C872 

17.531406403 

252.824924469 

-c . 

118. 740366936 

4.336146918 

11 5. 6 64155 96C 

11.337471962 

104.326683558 

-c . 

121.448160172 

4.4C 1327922 

65.653623581 

8.989222527 

56 .664 40 2 C 5 4 

-c. 

122.650173187 

4 . 423922424 

52.680656910 

8.311931133 

44.366725777 

-c. 

123 . 02 2 93 7 7 75 

4.466526926 

43.565114975 

7.807105064 

35.758009911 

-c . 

123. 316404343 

4.499121426 

36.918286324 

7.418309927 

29.499976357 

-c. 

123.553673744 

4. 52171 5920 

31. 923838654 

7.110941410 

24.812897444 

-c . 

123.749719620 

4.662C92927 

2C. 659891129 

6.344658494 

14.315232625 

-c . 

124.305338860 

4. 8 5766C949 

14.041173535 

5.810591578 

8.230582356 

-c. 

124. 769424438 

5. 1 636C 5969 

1C. 271 516747 

5.45149 1058 

4.820425685 

-c. 

125. 194414139 

5.6CC669046 

8. 51 5005775 

5.284209371 

3.234796405 

-c. 

125.649219513 

6. C 65477531 

8.354503302 

5.288650870 

3.106252422 

-c. 

125. 801103592 

6.6*4422466 

5.6C1751 924 

5.391392946 

4.210356577 

-c. 

126.229794502 

7. 129C4197C 

12.103486419 

5.425990999 

6.677455420 

-c . 

126.451245308 

7.22391 1464 

14. 13581C681 

5.430888414 

8.708922267 

-c. 

126.552917480 
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7.2713**211 

15. 733036160 

5.432405829 

10.300630221 

c • 

126.615120888 

7.2 1876C9S9 

18.000871 897 

5.434687614 

12.566184262 

-c. 

126.689459801 

7. 2 6 £2 1 51C6 

21.385763884 

5.440934181 

15 .9440 2570 2 

-c. 

126. 781545639 

7. 4 1 26 5 C4 53 

26.756425951 

5.459326029 

2 1 .2990955 2 2 

-0. 

126. 901064873 

7.514575461 

57. 415153C27 

5.701188564 

51.713964462 

-c. 

127.394231796 

7 • 5 5C9C 6446 

88. 815578461 

6.077634811 

82.7 37943645 

c • 

127. 718454361 

7. 561165527 

142.551 823576 

6.869661331 

135.722162241 

-c. 

128. 155666351 

7. 6C54C7S53 

229.236686578 

8.365262985 

22 0.8716229S2 

-c . 

128. 724721909 

7.623574436 

347.655506677 

10.682868950 

336 .973027720 

-c. 

129.390148163 

7.6*1740578 

538.311595731 

14.927673340 

523.383926 35 2 

-c . 

130. 414619446 

7.6*5014563 

64C.61750C3C5 

17.499565125 

623.117935161 

-c. 

131. 02 1 064 75 8 

7.661925572 

769. 2651556C7 

21.364532471 

747.900627136 

-c. 

132.105680466 

7.6 7CGGC017 

800.368080139 

23.048080444 

777 .319999655 

-c. 

132. 909090042 

7.676C74GC2 

77C. 538581848 

23.420906067 

747.117675761 

-c. 

133. 710794449 

7.6522C3462 

617. 803195768 

21.351158142 

596.452041626 

-c. 

134. 945844650 

7.654221573 

552.154174805 

20.899909973 

571.25426482 2 

c • 

135. 098976135 

7.7 1C37C004 

405. 576 62 1 899 

17.213455200 

388.763366695 

-c. 

136.104303360 

7. 720536467 

265.135163571 

14.023357391 

251.115806560 

-c. 

136. 856510162 

7.752756443 

161. 656569290 

11.431249619 

150.425319672 

-c. 

137.484615326 

7.763025534 

98.972272324 

9.691202164 

89 .282070 16C 

-c. 

137. 951662064 

7.015366556 

63. 102567739 

8.601566315 

54 .501401424 

-c. 

138.286458969 

7.621415567 

55. 73C446615 

8.363546371 

47.366900444 

-c. 

138.365261078 

7.655157310 

44. 955426178 

8.004860401 

36.990565777 

-c. 

138.492 776871 

7.670914654 

37.523616791 

7.744807005 

29.778005166 

c • 

138.593425751 

7.SC262 1956 

32.127161580 

7.550552368 

24.576605612 

-c. 

138.675113678 

7.9263*5342 

28.105384298 

7.402010202 

20.707374056 

-c. 

138. 742 950439 

7.5 5C066666 

25. 041463375 

7.286328793 

17.755134562 

-c. 

138. 800407410 

7.5737 64C25 

22. 646416C42 

7.194965124 

15.45345C5 16 

-c. 

138. 849868774 

7.5575C1373 

20.747635656 

7.122049689 

13.625585961 

-c . 

138.893045425 

6.C212 16856 

19. 21415SC12 

7.063432097 

12.150726514 

-c. 

138. 931209564 

8.1 16C6635C 

15.323706665 

6.922232032 

8.401474822 

-0. 

139. 052074432 

8. 3C5827375 

12.158560514 

6.857441485 

5.301115020 

c • 

139.210536957 

6. 76C174851 

1C. 456 8C8090 

7.092889428 

3 .3639 1862 2 

-c . 

139. 451656342 

8.565912840 

1C* 3 5969C151 

7.254678786 

3.145011336 

-C. 

139.521244049 

10.553345661 

12. 851 676578 

9.690242171 

3.201634316 

-c. 

140. 048994064 

12.322257931 

20.011977911 

15.569681168 

4.442296743 

-c. 

140. 624452591 

13.6 15011875 

31. 758782625 

25.375525236 

6.383257210 

-c. 

141.164415359 

14.6*7562365 

50. 9758541 C7 

41.632131577 

9 .347722520 

-c. 

141.735 769272 

15.455427332 

83. 721455389 

69.468955040 

14.252503753 

-c. 

142.401521683 

15.621124763 

112^295755151 

88.873520851 

23.426274300 

-c . 

142. 800090790 

15. 6 5C 527328 

143.081237793 

97.462361336 

45.618875560 

-c. 

142.930261612 

15.5C1812434 

146. 12152 7261 

98.715823174 

47.406104068 

-0. 

142. 962102 890 

15.524427354 

132. 05 706 9778 

98.347208023 

33.709861755 

-c. 

143. 045139313 

15.55C6245C5 

122.858665428 

100.063799858 

22.834865222 

-c. 

143.144659042 
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16. *<241*382 

1 78. 71 875C0C0 

151.763008118 

26.95574069C 

-C. 

143. 913 54 7516 

16.9742*9*Cl 

29C. 95941 5436 

249.589328766 

41.370065235 

-C. 

144. 896333694 

17.2*2995681 

425.926296234 

367.597076416 

58.329218266 

-c. 

145.797342300 

17.45C499735 

572.630378723 

496.161808014 

76.4685668C2 

-c. 

146.597604752 

17. 6 <25 831 13 

7 58, 4C65551 76 

659.242706299 

99.163848677 

-C. 

147.457208633 

17.7*5833139 

914.821815491 

796.712173462 

118.1096363C7 

-c. 

148. 090398788 

17*6 *9G 631*6 

1122.873321533 

979.727584839 

143.145726*54 

-C. 

148. 847354889 

17.92791*517 

1301. 967193604 

1137.384201050 

164.582968725 

-C. 

149.438930511 

ie.CC*7498*8 

1 52 5. 7C2667236 

1334.453704834 

191 .24895668C 

- C . 

150.120365143 

18.C 7556322C 

leOS. 793594360 

1584.829635620 

224.963947256 

-c. 

150. 914340973 

18. 12*999645 

2132. 767639160 

1869.627593994 

263.1400451*6 

-c . 

151.742181778 

18.19C9997A* 

2490.279571 533 

2185.030517578 

305.249027252 

-c. 

152.587135315 

18.217999697 

27C2. 114196777 

2371.975585937 

330.138595561 

-c. 

153. 058361053 

18.2<499<*47 

2940. 828887939 

2582.689117432 

358.139751434 

-c. 

153.568162918 

18.271999836 

321C. 970153809 

2821.195739746 

389.7743873*0 

-c. 

154.121232986 

16.298999766 

3516.C2941 8945 

3092.357147217 

425 *672260 264 

-c. 

154. 723484039 

18.325999737 

3868.69046C205 

3402.091583252 

466 .598850250 

-c. 

155.381383896 

18.352999667 

4271.116577140 

3757.628692627 

513.487876852 

-c. 

156. 103029251 

16.379999638 

4735.325622559 

4167.841796875 

567.4837795C7 

-c. 

156. 898090363 

16.4C6999626 

5273.666625977 

4643.671203613 

629.9954223*2 

-c . 

157.777486801 

18.43399977? 

59C1.38885498C 

5198.629089355 

702.7597503*6 

-c. 

1 58. 754465103 

18.4 6C999727 

*637.367878418 

5849.462524414 

787.92532111* 

-c. 

159. 845474243 

18.4*7999676 

75C5. C30822754 

6616.886474609 

888.144305556 

-c. 

161. 070623398 

18. 514999626 

8533.066C59570 

7526.411987305 

1006.674041748 

-c . 

162.453268051 

18. 5*199981 7 

9756. 593C17578 

8609.121826172 

1147.471176147 

-c. 

164. 022 789001 

18.5*69997*7 

11217. 256956008 

9902.029052734 

1315.227890015 

-c. 

165. 815568924 

16. 595999718 

12962.794799805 

11447.509643555 

1515.285095215 

-c . 

167. 872980118 

16. *229996*8 

15C43. 536376953 

13290.308349609 

1753.227920522 

-c . 

170.246139526 

18. **9999857 

17503.840087891 

15469.953125000 

2033.886871238 

-c. 

172. 990770340 

16.656999661 

18412.998779297 

16275.582031250 

2137.41665*454 

-c. 

173. 998075485 

18.6(7999744 

19366. 775634766 

17120.853271484 

2245.922149*58 

-c. 

175. 056886673 

18.6769998C7 

203*4.027(99*09 

18004.770019531 

2359.257080076 

-c . 

176.170106888 

18.6 659996 7C 

21402.589843750 

18925.431640625 

2477.158172*07 

-c. 

177.336496353 

18.65*999*95 

22479. 0C7324219 

19879.793701172 

2599.213592525 

-c. 

178.562240601 

18. 7 C3999756 

23588.414794 922 

20863. 564453125 

2724.850219727 

-c. 

179. 844430923 

18.721999864 

25877.522705078 

22893.997070312 

2983.525572720 

-c . 

l 82.591482162 

18. 73C9997C6 

27037.84887*953 

23923.527587891 

3114.321075425 

-c . 

1 84. 059986114 

16. .739999771 

28192.196777344 

24948.007568359 

3244.189208564 

-c . 

185.587839127 

18. 7*6999834 

29325.350585937 

25953.965820312 

3371.3847351C7 

-c. 

187.178663254 

16.757999897 

30420. 053955C78 

26926.118896484 

3493.934967C41 

-c. 

188. 826787949 

16.7**999722 

31457.287353516 

27847.615234375 

3609.672058105 

-c. 

190. 533632278 

16.775999764 

32416.898193359 

28700.591064453 

3716.306945601 

-c. 

192.290845871 

18.7 64999847 

33276.187500000 

29466.691650391 

3811.495727525 

-c. 

194. 097133636 
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18. <«29S991C 

34C2C. 933105469 

30127.974853516 

3892.95816C4QC 

-c. 

195 . 94 3 47 1 90 9 

18.8C2995725 

34626. 4C2343750 

30667.810302734 

3958.59185791C 

-C. 

197. 023368073 

]8.eiiscc 7 c £ 

35078. 497558594 

31071.894042969 

4006 .60 3 3020 C 2 

-C . 

199. 73 0163 5 74 

18.8209996*1 

35364. 700683554 

31329.087158203 

40 35.613432836 

-c. 

201.653753281 

18.829999924 

35476. 984663281 

31432.224609375 

40 44 . 7 600 70 GC 1 

-c. 

203. 585494995 

18.828999746 

35412.376464844 

31378.617431641 

4033.758819580 

- C • 

205.515291214 

18.8*7999811 

35173.141601562 

31170.218750000 

4002.922668457 

“ c • 

207.434228897 

18. 66*999674 

34766. 63623C465 

30813.495117187 

3953.141021729 

-c. 

209.331714630 

18. 8*5999699 

34204. 76 5042 96 9 

30318.951416016 

3885 .8175354CC 

-c . 

211.200071335 

18. 6749997*2 

335C3. 124511 719 

29700.357910156 

3802.766540527 

-c. 

213.032451630 

18.6 62999625 

3268C. 024414062 

28973.917480469 

3706.106923594 

-c. 

214. 821762085 

16.692999*49 

31755. 42C41C156 

28157.286132812 

3598. 134155272 

- c ■ 

216. 561702728 

18. 9C 19997 12 

20749. 756347656 

27268.565917969 

3481.1902160*4 

-c . 

218.245018005 

18.91C999775 

29663.181 684766 

26325. 605468750 

3357.576385456 

-c. 

219.871004105 

16.919999599 

26574. 738C371C9 

25345.27758789 1 

3229.460296*21 

-c . 

221.436441422 

18. 9*699 9766 

25161.213378506 

22324.641845703 

2836.57144165C 

~~ c • 

225. 761177063 

18. 9 5 59996 12 

24037. 8e330C78l 

21330.169921875 

2707 .7 13 1956G 1 

-c. 

227. 074735641 

18.9*499967* 

2293e.251Cl 5625 

20356.536132812 

2581.754760742 

-c. 

228.328451157 

16.972999729 

21869.373535156 

19409. 906494141 

2459.46697556C 

-c. 

229. 523195267 

16.9629995*2 

2C836. 324951172 

18494.901855469 

234 1.422882G6C 

-c . 

230.660207748 

18.991999626 

19642. 7517C8984 

17614.736572266 

2228 .0 1498*131 

-c. 

231. 742801666 

19.CCC9 99689 

18851. C642C6564 

16771.562988281 

2119.501096633 

-c . 

232. 772596359 

19.CC 99 99752 

17562. 551552734 

15966.575073242 

2016.016432716 

-c . 

233. 748632431 

19.0269 997C2 

15517. 536885648 

13782. 168457031 

1735.76832 56C6 

-c . 

236.408481598 

19. C <29996 53 

13426. 1C852C5C8 

11927.599575586 

1498.508499 146 

-c. 

238.698686600 

19.C 9C9996C3 

11665.651577539 

10366.358886719 

1299 .29301452* 

-c. 

240.676637650 

19. 1 17999554 

1C186. 147216757 

9055.687377930 

1132.459777832 

-c . 

242.394172668 

19. 1**9995C4 

8547.1 C 82 7636 7 

7954.488098145 

992.6201625*4 

-c. 

243. 893074036 

19.171999693 

7501. 553C35551 

7026.504516602 

875 .0 4848 4602 

-c. 

245.208374023 

19.1989996*3 

7 Cl 6. 915128418 

6241 . 145080566 

775.774022553 

-c. 

246. 369470596 

19.255999594 

6264.680175781 

5573. 151672363 

691.528457*42 

-c. 

247. 399978638 

19.2^2999544 

5621. 553566695 

5001.938293457 

619.6552658C8 

-c. 

248.319984436 

19.279999495 

5C68. 815734663 

4510.8 10668945 

558.005025400 

-C. 

249. 144590378 

19.26*999663 

4591. C648388b7 

4086.244995117 

504.83983612 1 

-c. 

249. 887767792 

19.322999624 

41 76.C38C24902 

3717.283905029 

458.754119873 

-0. 

250. 560745239 

19.2 <C9 9 95 64 

3613.61C778809 

3395.005401611 

418.605361938 

-c . 

251.172679901 

19. 267999535 

3495.589172363 

3112.126525879 

383.462642670 

-c . 

251. 731212616 

IS .4 149 9 9465 

3215.239471436 

2862.675842285 

352. 5636 1770* 

-c . 

252.243165970 

19. 4*1999674 

2967. 01 C58S6C0 

2641.729736328 

325.280841827 

-c . 

252. 714048386 

19.4*6999624 

2746. 3121643C7 

2445.216033936 

301.096130371 

-c. 

253. 148424149 

19.495999575 

2545.317718506 

2269.739501953 

279.578216552 

- c • 

253.550615311 

19.5*5959714 

2214.158721924 

1970.997711182 

243.161003112 

-c. 

254.273469925 

15 .6 1*465622 

1865.67C6S5306 

1678.049087524 

207.8218076*1 

-c. 

255. 038740158 

19.662925525 

1 626.685791 C16 

1446.402618408 

180 .283166665 

-c . 

255.695110321 
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19.771566153 

1356.706535697 

1206.323348999 

152.383584576 

-c. 

256. 442337036 

15. 8*2349^1 

1111.083523340 

983.485603333 

127.598315239 

-c. 

2 57.224601746 

2C.0152e<54* 

895.655387878 

791.468493230 

108 • 18688564 £ 

-c. 

258. 010299683 

2C. 170386314 

73C. 765558840 

634.734779358 

96.031217575 

-0. 

258. 798522949 

20.2*7635561 

6C8.25844574C 

513.500259399 

94.758184423 

-C. 

259.633174896 

20.524053045 

550. 6C361 4 607 

437.741840363 

112.861772537 

-C. 

260.533470154 

20 • 56S2C 6238 

540. 6C5878540 

425.843894958 

122.765981674 

-c. 

260. 787609100 

20 .615585576 

575.936172485 

413.169872284 

166.766298254 

-c. 

261.565162659 

20.746455723 

642.285335355 

424.291118622 

217.994220734 

-c. 

262.182483673 

2C • 7507 7 59 15 

715. 200851440 

446.215858459 

272.984952561 

-c. 

262.706447601 

20. 8 250 * 6346 

847. 166603C88 

488.789493561 

358.377109526 

-c. 

263.378993988 

50.857243061 

542.172332764 

522.575355530 

419.596977234 

-c. 

2 63. 793 087006 

20. 667257153 

1124. 365384766 

589.885032654 

534.484352112 

-c . 

264.479995728 

20.5 1C5 65761 

1334.332092285 

669.707359314 

664.624732571 

-c. 

265.159774780 

20.556764754 

1523.617156582 

742.859634399 

780.757522563 

-c. 

265. 706092834 

50.5*2575746 

1765.64122C093 

839.007362366 

930.63385772? 

-c. 

266.351768494 

50.550477153 

1520. 679531641 

898.462272644 

1022.21765855? 

-c. 

2 66. 72 1 8093 8 7 

50.556374735 

2055.258392334 

967.492752075 

1127.765640255 

-c. 

267. 130096436 

20.566272116 

2258.064880371 

1048.022476196 

1250.042404175 

-c . 

267.580207825 

2C. 5 7416572 1 

2534. 92 C 667920 

1142.444717407 

1392 .476 15C 5 1 3 

-c. 

268.081710815 

20.562C671C6 

2813.0C357C557 

1253.716217041 

1559.287353516 

-c. 

268. 640037537 

20.565564724 

3141.227294522 

1385.517410278 

1755.709884644 

-c. 

269.268959045 

20.5578621C1 

3530.525509424 

1542.375396729 

1988.150512655 

-c. 

269. 976799011 

21.0C5755716 

3S54. 261657715 

1729.840393066 

2264.421264646 

-c. 

270. 78363 0371 

21.C12657C53 

4548.258522549 

1954.530426025 

2593.768056524 

-c. 

271. 702091217 

21.0215547C6 

5210. 842346191 

2224.085876465 

2986.756465727 

-c. 

272.761329651 

21.0 5545 5C 65 

6001. C3009C332 

2546.617614746 

3454.412475566 

-C . 

273. 979446411 

21.0273457C1 

6935. 812866211 

2929.461639404 

4006.351226807 

-c. 

275.396488190 

21.0*2614656 

7643. £55468750 

3220.310607910 

4423 .54486 C8 4 C 

-c. 

2 76.457962036 

21.C47675656 

8417.766235352 

3539.054870605 

4878.711364746 

-C . 

277. 623435974 

21. C! 2144653 

5248. 875CCOOOO 

3882.361083984 

5366.51391601 6 

“ c • 

278.929851532 

21 .C 626746 6 8 

10995. 751553125 

4609.543945312 

6390.208007812 

-c. 

281. 879650116 

21.066925666 

11847.1 75503320 

4964.089233398 

6883.086665522 

-c. 

283. 546550751 

21.0715723C4 

12241. 763305664 

5130.050720215 

7111.712585445 

-c. 

284. 432365417 

21.C742C4662 

12606.785278320 

5284.289306641 

7322.495971680 

-c. 

285.349369049 

2 l.C 16 8 273C 1 

12534.436523437 

5423.573974609 

7510 .862546826 

-c. 

286. 283180237 

21 .C 75465681 

13216. 556665C39 

5544.675292969 

76 72 . 28 1 37 207C 

-c. 

287.241283417 

21.C621C2259 

13447.266113281 

5644.631774902 

7802.634336375 

-c. 

288.211193085 

21. C (47 246 76 

13615.175321289 

5720.833374023 

7898.345547266 

-c. 

289.196231842 

21.067367256 

13727. 562402344 

5771.257446289 

7956 *7 04956 C 5 5 

-c . 

2 90.186969757 

2 1. C 655 5 5514 

1377C.5357666C2 

5794.551940918 

7975.983825664 

-c. 

291.181663513 

21.052632254 

13745.764160156 

5790.159667969 

7955.604492187 

~c. 

292.175117493 

21.055264512 

13654. 5C2C75195 

5758.342407227 

7896.159667565 

-C . 

293.164260864 

2 1 .057857251 

13499. 56CC58594 

5700.169250488 

7799 .39080 8 1C 5 

-c. 

294. 141128540 
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21.1CC52S9C9 

13285.433471680 

5617.412109375 

7668 .C 2 1 36 2 20 5 

-c . 

295. 091625214 

21. 1C rl622€«= 

13C18.C992431 64 

5512.461303711 

7505.637939452 

-0. 

296. 035865784 

21. 1CS7S4SC7 

127C4. 477905273 

5380.111206055 

7316.366699219 

-c. 

296. 941566467 

21. 1(8427266 

12352. 2C8374C23 

5247.463256636 

7104.745117 lei 

-c . 

297. 840232849 

2I.1]1C‘{<SC4 

11565.027343750 

5093.676452637 

6875.350851 1 12 

-c . 

298. 691669464 

21. 1 H324SC2 

11140.267475586 

4759.136901855 

6381.150572720 

-c. 

300. 325843811 

21. 15^22227? 

9839. 534301 758 

4230.496154785 

5609.038146572 

-c. 

302.558078766 

21. 129467276 

8992.31 5185547 

3884.463012695 

5107.852172852 

-c. 

303. 887386322 

2 1 • 1 34752274 

8151.074462891 

3556.170227051 

46 3 4 .904 23 5 8 40 

-c . 

305. 100105286 

21. 1*C0 17 27 1 

7445. 44 e6 6 5434 

3251.466491699 

4197.982177724 

-c . 

306.177421570 

21.142645659 

71C3. C29113770 

3108.867279053 

3994.161824717 

-c . 

306. 681 842804 

21.1 505 <7266 

6163.052795410 

2721.064056396 

3441 .9887290 14 

-c . 

308. C6 91 33 759 

21. 1564<4661 

5364. 54561 7676 

2390.738433838 

2974.207182826 

-c . 

309.258323669 

21. 166342258 

4652. 9981C7510 

2111.798889160 

2 5 81. 199 2 1 8 7 5 C 

-c. 

310.289493561 

21.174225874 

4128. 7 5C3051 76 

1876.899139404 

2251.851165771 

-c. 

311.181037903 

21.162127251 

3654. 54858C713 

1678.936935425 

1975.612045268 

- c . 

311. 561 956024 

21. 1 SC 0 34 6 66 

3254. 771728516 

1511.585174561 

1743.186552555 

-c. 

312. 654319763 

21. 157922243 

2516.281707764 

1369. 500061035 

1546.78 1646725 

-c. 

313.267261505 

21.2(5825655 

2628.2C2 72 8271 

1248.241577148 

1379 .96115 1 122 

-c. 

313. 809421539 

2 1.2 12727226 

2381. 7C4162558 

1144.195236206 

1237.508926252 

-c . 

314.294589996 

21.221624851 

21 65. 557686768 

1054.41162 1094 

1115.186065674 

-c . 

314. 728404999 

21.2 <5522226 

1586.086257535 

976.507141113 

1009.579116821 

-c. 

315. 120410919 

21. 227415844 

1826. 44812C117 

908.536895752 

917.911224265 

-c. 

315. 478000641 

21.252214636 

1564.151986084 

796.373573303 

767.818412781 

-c. 

316.104190826 

21.265CC5625 

1355. 8C7C83130 

708.403335571 

651.403747555 

-c. 

316. 631237030 

2 1. 2527C2156 

1128. 545523526 

608 .212631226 

520 • 7368527CC 

-c . 

317.2 83466339 

21. 3CC5<5813 

1C67.0C859C698 

581.132408142 

485.876102556 

-c. 

31 7.469886780 

21. 4 2^3 2 7 1 2 5 

565. 0227127C8 

355.546531677 

209 .476 175123 

-c. 

319. 483890533 

21. 5*8054655 

39C. 16C3C8e38 

270.669521332 

119.490786552 

-c . 

320.4312705 99 

21.671762255 

3C4. 751793823 

225. 586166382 

79.205626466 

-c. 

321.000110626 

21.5 15237127 

220. 1 81337357 

175.728292465 

44.453044415 

-c . 

321.702465057 

22.41*147125 

148. 271181107 

126.238922119 

22.032256567 

-c. 

322.445617676 

23 .26C22582C 

55. 953726730 

84.941649437 

1 1.052076657 

-c . 

323. 074085236 

24.46C35717C 

64.471681595 

57.667510033 

6.804 171 1*5 

-c . 

323. 523757935 

25. 721375747 

51. 966CO1034 

43.810389042 

8.155611752 

-c. 

323. 896274567 

26.267365755 

55.483874321 

43.058216572 

16.425657745 

-c . 

324. 158714294 

26.42635 2C78 

67. 555 567 7 C3 

44.999886513 

22.595681150 

-c. 

324.261611938 

26. 56535464C 

84. 3C485 1 532 

49.878439903 

34.426411625 

-c. 

324.411190033 

26.7C43572C2 

124.16 5C53 368 

62.772307873 

61.392745455 

-c. 

324.661231995 

26.8 52272C34 

266.025745352 

112.275869370 

155 .749876022 

-c. 

325.260646820 

26.656624575 

383. 95513C249 

153.137489319 

230.861640520 

-c . 

325.579940796 

26.5232 521C6 

458. C57495117 

193.672779083 

304.384716024 

-c . 

325. 845169067 

26.545865633 

673.436627515 

256.400641815 

417.0 30185655 

-c . 

326. 202972412 

26.5676 14651 

645.919326782 

318.432830811 

527.486495572 

-0. 

326. 515129089 
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26.9164E716C 

956.471397400 

358.315246582 

598.156 15C616 

-C. 

326.700984955 

26.965359669 

1C88. 019134521 

405.879386902 

682.139747650 

- C . 

326. 915142059 

26.99*222176 

1245.037811279 

462.786026001 

782.251785576 

-C. 

327. 159030914 

27.0C31C4667 

1432. 640197754 

530.943305969 

901.696891785 

-c. 

327.438171387 

27. C 11977196 

1656.214294434 

612.381286621 

1043.833007815 

-C. 

327.759063721 

57. C238C7C49 

2018.130401611 

744.658325195 

1273.472076416 

-c. 

328.268108368 

2 7.C2562714C 

2454.659454346 

904.925064087 

1549.734390259 

-c. 

328. 894817352 

27.0*0424576 

3067.539001465 

1131.483871460 

1936.055130005 

-c. 

329. 857528687 

27. C £9297323 

3417. 4CS17S687 

1262.056610107 

2155.352569560 

-c. 

330.530647278 

27.C6S2 1225C 

3612.253204346 

1335.568542480 

2276.684661865 

-c. 

331.018394470 

27. C 711271 76 

3756.208129883 

1390.810760498 

2365.397369365 

-c. 

331.526832581 

27.077042341 

3834. 566589355 

1422.288970947 

2412.277616408 

-0. 

332.050621033 

27.C79999924 

3846.136932373 

1427.967773437 

2418.169156936 

-c. 

332.314563751 

27.CE591485C 

3812. 518554687 

1418.533325195 

2393.985259495 

-C. 

332.839714050 

27.C9182977? 

3707.567 Cl 66 02 

1382.829040527 

2324.737976C74 

-c. 

333.350139618 

27.CS77447C3 

3542.947021484 

1324.987396240 

2217.959625244 

-c. 

333. 840145111 

27.1C6617212 

3220.007232666 

1209.670166016 

20 10 .3370666 5C 

-c. 

334.53 0418396 

‘27. 1 164473C4 

2733.355C72021 

1033.961044312 

1699.394027710 

-C. 

335. 337635040 

27. 12C27739S 

2268.305938721 

864.81050 1099 

1403.495437622 

-c. 

336.010002136 

27. 129149SC4 

1961.611297607 

752.742416382 

1208.868881226 

c • 

336.429584503 

27. 1*6022175 

1696. 824554443 

655.678359985 

1041.146194456 

-c. 

336. 796688080 

27. lit e«46E4 

1472.147262573 

573.086502075 

899.060760496 

-C . 

337.110435486 

27. 165767193 

1283.1C55SC820 

503.417854309 

779.687736511 

-c. 

337.382797241 

27.1 746397C2 

1124.478973389 

444.820327759 

679.658645630 

c • 

337.620277405 

27. 1E3512211 

991.261291504 

395.499275208 

595.762016296 

-c. 

337. 82 83 843 99 

27 • 2C 12 f 7 229 

784.1C623S319 

318.567932 129 

465.538307190 

-c. 

338.172492981 

27.245619774 

479.259326935 

204.623519897 

274.635 80 7C37 

-c. 

338.759632111 

27. 26 11C981C 

34 9* 4 75C67139 

155.680753708 

193.794313431 

-c. 

339. 063842773 

27.3 16599846 

268.565544128 

124.934301376 

143.631242752 

-c. 

339.2 83172607 

27.4CC769711 

166.891721725 

85.819760323 

81.071961403 

-c. 

339. 62 8677368 

27.464929814 

119.265119553 

67.120431900 

52.144687653 

-c. 

339. 832988739 

27. 5 6 9 1 C 96 76 

93. 153319359 

56.635271549 

36.518047810 

-c. 

339. 968540192 

27.737449646 

66. 732913017 

45.641374588 

21.091538429 

-c . 

340. 143898010 

27.969959717 

49.964362144 

38.152446747 

11.811915159 

-C. 

340.293037415 

28.579149723 

36.361289024 

31.114526510 

5.246762335 

-c . 

340.47094 7266 

30.CCC0CC0QC 

26. 717815399 

24.564364195 

2.15345 12C4 

-C. 

340.651081085 


FOR 115 RESONANCES, 345 ENERGY PCINTS ARE REQUIRED FOR DEVIATIONS OF LESS THAN 
10.0 BARNS AND 5.0 PERCENT AT INTERPOLATED POINTS. 

OF THE INITIAL 861 PCINTS, 548 HERE OISCARCEC AND 32 ACCIT10NAL POINTS WERE INSERTED. 
IS POINTS WERE RETAINEO AS WAX, PIN CR INFLECTION. 


AVERAGE PERCENT DEVIATION CF DISCARDED PCINTS 

1.97583 
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